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ABSTRACT 


Simultaneous  measurement  of  the  extinction  coefficient  ae  and  liquid  water  content 
W  were  made  for  laboratory  generated  cloud  at  C02  laser  wavelength  X  *=  10.591  micro¬ 
metres  (pm).  Good  agreement  between  measurement  and  prediction  (Chylek,  1978)  was 
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found  for  a  wide  range  of  W  (up  to  7  g  m  )  and  oQ  (up  to  1  m  ). 

The  first  definitive  measurements  of  backscatter,ob>  and  extinction  coefficients 
made  on  laboratory  -  generated  fog  droplet  distributions  at  visible  (X  =  0.6328  pm) 
wavelengths  were  in  good  agreement  with  a  size  distribution  independent  relation 
j0  =  17.7  ab  (Pinnick,  Jennings  et  al,  1983). 

The  first  definitive  measurements  of  backscatter  coefficient  and  extinction  coeffic¬ 
ient  made  on  laboratory  cloud  are  reported  at  C02  Laser  wavelengths.  The  measure¬ 
ments  yieldinq  o^/o ^  ratios  of  between  350  and  about  550  are  in  good  agreement  with 
numerical  calculations  performed  on  narrow  size  distributions  of  natural  cloud. 
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STATEMENT  OF  RESEARCH  PROBLEM 


Research  conducted  during  the  term  of  this  contract  consisted  oc  measurement 
of  transmission  and  backscatter  of  electromagnetic  (em)  energy  through  water 
droplet  clouds  with  k'-own  size  distribution  and  liquid  water  content  at  visible 
and  infra-red  (IR)  wavelengths.  Theoretical  analyses  o+  relations  between 
extinction,  absorption  and  backscatter  which  incorporate  size  distribution  de¬ 
pendencies  were  also  made. 


SUMMARY  OF  MOST  IMPORTANT  RESULTS 

Experimental  measurements  made  of  extinction  coefficient  o  ,  m’"*  and  liquid 
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water  content  W,  g  m  for  laboratory  aenerated  cloud  at  the  CCt,  laser  wave¬ 
length  A  =  10.591  micrometres  (urn)  verified  the  size  distribution  independent 
relation  ae/W  =  3ttc/2A  derived  by  Chylek  (1978).  The  first  continuous  set  of 
measurements  of  c  and  W  are  reported. 

The  first  definitive  simultaneous  set  of  measurements  of  volume  backscatter  and 
extinct' on  coefficients  o^o  for  laboratory  cloud  at  visible  wavelengths 
(A  =  0.6328  pm)  are  reported.  The  measurements  show  qood  agreement  with  the 
theoretically  predicted  size  distribution  independent  relation: 

Ge//°b=  8~/g(X),  (Pinnick  et  al ,  1983),  where  g  (A)  is  a  slowly  varvina  function 
of  wavelencth. 

The  first  definitive  set  of  measurements  of  volume  backscatter  coefficient  a  and 

o 

volume  extinction  coefficient  ae  for  laboratory  generated  cloud  at  CO2  Laser 
wavelengths  (10.261  -  10.591  pm)  are  reported.  The  measurements  of  extinction 
to  backscatter  ratio  ^e/o^  possessing  values  ranging  from  350  upward  give 
relatively  aood  agreement  with  numerical  calculations  carried  out  on  relatively 
narrow  size  distributions  for  cumulus  and  stratus  cloud  (Pinnick  et  al,  1983). 
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Electro-Optical  Transmission  and  Liquid  Water  Content 
of  Fogs  and  Clouds. 


Introduction: 

A  possible  relation  between  extinction  and  liquid  water  content  is  of 
considerable  interest  and  has  preoccupied  atmospheric  scientists  over  the 
past  several  decades.  The  determination  of  simple  relationships  between 
electro-optical  (E-0)  transmission  and  measureable  meteorological  parameters 
such  as  cloud  water  content  (LWC)  could  be  used  to  predict  transmission  at 
middle  infrared  (IR)  wavelengths  from  a  knowledge  of  cloud  or  fog  LWC. 
Conversely,  a  measure  of  the  transmission  at  middle  IR  wavelengths  could 
yield  a  value  of  liquid  water  content  using  validated  relations  from  this 
work . 

Most  previous  work  has  been  directed  towards  seeking  relationships  between 
liquid  water  content  and  visibi 1 i ty,  for  example: 

Eldrlcige  (  1966,  1571  ),  Barteneva  and  Polyakova  (1965)  and  Kumai  (1973). 
Pinnick  et  al  (1978)  give  approximate  empirical  relationships  between 
extinction  coefficient  .o  and  liquid  water  content  VJ  based  on  the  work  of 
the  forementioned  authors.  However,  Pinnick  et  al  (1978)  found  that  the 
predicted  values  of  visible  extinction  o  differed  by  about  an  order  of 
magnitude  for  the  same  LWC  (0.01  g  m  ).  Thus  any  size-distribution  - 
independent  relation  between  extinction  at  visible  wavelengths  and  LWC  cannot 
be  applied  to  fogs  in  general. 

Extinction,  Absorption  and  Liquid  Water  Content  of  Water  Clouds 

Consider  a  pol yd  is  pension  of  spherical  water  droplets  described  by  the  size 

distribution  n(r),  where  r  is  the  radius  or  a  given  droplet.  We  examine 

relationships  between  the  extinction,  r  ,  absorption  coefficients ,  a  ,  and 

e  a 

the  liquid  water  content  W  given  by 

°e  =  /r2Qe(n,.x)n(r}dr  (!.'} 

'e  '  J  ’L 0a (n, x )n{ r )dr  (1.2' 

--  3 dr  (1.3) 

wnere  ,  the  liquid  droplet  density,  Qe  (m,x) ,  Qa  (m,x)  are  the  efficiency 
factors  for  extinction  anc  absorption  for  a  droplet  with  refractive  index  m 
and  size  parameter  x,  defined  by  the  ratio  of  the  droplet  circumference  to 
radiation  wavelength  .  It  has  been  shown  by  Chylek  (  1978)  that  the 
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efficiency  factor  Qe  (m,x)  can  be  approximated  by  a  linear  function  of  droplet 
size  parameter 

Qe  (m,x )  =  c  (X)  x  (1,4) 


providing  size  parameter  x  l  xm  (xm  =  2TTrm/A).  Pinnick  et  al  (1979)  have 
also  shown  that 


Qa  (m,x)  =  c 1  (A)  x  (1.5) 

The  use  of  these  simple  linear  relationships  for  the  Mie  efficiency  factors 
in  the  expressions  for  the  extinction  and  absorption  coefficients  given  by 
Eqs.(l)  and  (2)  reduce  these  coefficients  to  the  simple  form: 


o 


e 


(1.6) 


a 


e 


(1.7) 


where  c  (X)  and  c1  (X)  are  the  slopes  of  the  straight  lines  approximating 
the  Mie  efficiency  curves. 


The  wavelength  at  which  these  linear  relationships  are  valid  is  determined 
by  the  radius  of  the  largest  droplets  present  in  the  polydispersion. 

Chylek  (1978)  and  Pinnick  et  al  (1979)  give  values  of  the  maximum  radius  rm 
and  corresponding  wavelengthsX,  required  for  the  validity  of  (6)  and  (7). 

They  find  that  the  relation  (6)  should  work  best  for  water  droplets  in  fog 
or  cloud  at  or  near  X  =  11  pm  where  rm  is  found  to  be  14  pm.  However, 
relation  (6)  might  also  be  expected  to  work  reasonably  well  for  fog  or  cloud 
droplets  for  the  wavelength  range  9.5 <  X  <  11  pm  since  12.5 <  rm <  14  pm 
at  these  wavelength  values.  One  might  expect  relation  (6)  to  work  reasonably 
well  even  if  droplets  in  a  particular  distribution  have  radii  greater  than 
14  pm  providing  they  do  not  contribute  excessively  to  either  the  liquid  water 
content  or  the  extinction. 


Pinnick  et  al.  (1979)  verified  theoretical ly  relation  (6)  within  a  factor  2 
at  A  =  11  pm  for  341  droplet  size  distribution  measurements  of  atmospheric 
fogs  formed  under  a  variety  of  meterologi cal  conditions,  and  for  which  reliable 
fog  data  was  available.  Pinnick  et  al.  (1979)  used  the  raw  size  distribution 
fog  data  of  Pinnick  et  al.,  (1978),  Garland  (1971),  Garland  et  al.,  (1973), 

Roach  et  al.  (1976),  Kunkel  (1971),  and  Kumai  (1973)  to  calculate  the  liquid 
water  content  and  the  extinction  coefficient  for  all  of  the  fogs,  and  compared 
the  results  to  the  theoretical  relation  (6). 

One  of  the  objectives  of  the  work  was  to  investigate  experimentally  the  validity 
of  relation  (6)  i.e,  oQ  -  .  W  by  making  simultaneous  measurements  of 


-  3  - 


(a)  the  extinction  coefficient  ae 

(b)  the  liquid  water  content  W 

and  (c)  the  size  distribution  of  the  droplets. 

This  approach  is  distinctly  different  from  that  of  Pinnick  et  al  (1979)  in 
that  the  extinction  coefficient  and  LWC  is  measured,  rather  than  calculated 
from  the  drop  size  distributions.  One  reason  for  measuring  the  droplet  size 
distribution  is  that  the  validity  of  relation  (6)  requires  droplets  must  be 
less  than  a  certain  size  as  discussed  earlier  in  1.2. 

1.3.  Experimental  Apparatus  and  Techniques 
3 

A  1  m  chamber  of  path  length  L  =  lm,  constructed  of  waterproof  plywood  was 
used  for  all  measurements.  The  chamber  walls  were  painted  matt  black  to 
minimize  stray  light  reflections.  Saturation  of  the  chamber  was  maintained 
through  the  use  of  matt  black  absorbing  cloth  material  lining  the  chamber  walls. 
A  5  cm  diameter  opening  at  the  front  side  of  the  chamber  permitted  the 
entrance  of  laser  radiation  (either  a  He-Ne  or  Co£  laser  beam)  into  the 
chamber  and  also  allowed  for  the  extraction  of  backscattered  light  from  the 
cloud  droplets  in  the  chamber.  The  laser  beam  traversed  the  chamber  path 
length  and  emerged  through  a  small  aperture  (  few  mm  diameter)  in  the  rear 
exit  wall  of  the  chamber 

Figure  1.1.  is  a  schematic  diagram  of  the  optical  arrangement  used  in  the 
extinction  coefficient  measurements.  It  consists  of  a  tunable  vertically 
polarised  CC^  laser  source  (Sylvania  Model  941 S )  with  wavelengths  available 
over  the  wavelength  range  10.2  to  10.7  micrometres  (pm).  The  CO^  laser 
transmissions  path  is  made  coincident  with  a  visible  (0.6328  pm)  He-Ne  laser 
transmission  path.  Alignment  is  accomplished  with  aluminized  mirrors  which 
are  mounted  on  micrometer-controlled  translational  and  rotational  stages. 

A  portion  of  the  CO^  laser  beam  is  reflected  by  a  ZnSe  window  for  a  laser 
reference  signal.  The  main  beam  enters  the  fog/cloud  chamber  of  pathlength 
L  through  a  narrow  window  and  exits  through  a  second  narrow  window  into  a 
second  dectector.  The  main  beam  and  reference  detectors  (Laser  Precision 
Corporation  Models  RkP-545  and  Rk-5100)  consist  of  a  pyroelectric  laser  probe 
together  with  a  synchronous  radiometer  readout.  The  synchronous  ratiometer 

readout,  model  Rk -5200  is  used  for  the  output  of  the  ratio  of  the  main  beam 

”6  ~2 

and  reference  signals.  The  detectors  have  negligible  drift  (about  10  W  cm 
after  warm-up)  compared  to  typical  radiance  levels  monitored  by  the  detectors 
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(of  the  order  of  1  W  cm  ).  The  radiation  wavelength  is  monitored  precisely 
with  a  spectrum  analyzer  (Optical  Engineering  Inc.  Model  16-A  CO^  Laser 
Spectrum  Analyzer). 

Appropriate  considerati on  is  given  to  optimum  transmissometer  design  to 
reduce  forward  scattering  corrections  -  in  accordance  with  the  results  of 
Deepak  and  Box  papers  a,  b,  (1978).  For  example,  the  use  of  a  10  cm  lens 
which  focuses  the  CO^  beam  of  half  width  2.5  mm  through  an  aperture  of  radius 
0.5  rnm  will  result  in  a  correction  of  less  than  2%  due  to  forward  scattering 
at  X  =  10.6  pm  for  a  Deirmendjian  haze  M  (broad  distribution)  of  mean  droplet 
radius  5  pm.  In  this  work,  aperture  diameters  from  1  to  2.5  mm  were  used. 

Water  droplet  clouds  were  generated  within  the  chamber  by  a  pair  of  commercially 
available  "cool-mist  vapouri zers" .  The  cloud  droplet  size  distribution  was 
determined  by  a  Particle  Measuring  Systems  (PMS)  classical  scattering  aerosol 
spectrometer  probe  (CSASP)  which  can  sense  water  droplets  with  radii  from  0.23 
to  14  pm  (Pinnick  and  Auvermann,  1979).  Corrections  were  applied  to  the 
measured  droplet  ,ize  distributions  in  accordance  with  the  response  calculations 
for  the  Knollenberg  light-scattering  CSASP  counter  (Pinnick  and  Auvermann,  1979). 
The  cloud  generator  covered  the  range  of  droplet  sizes  usually  encountered  in 
natural  cloud.  The  cloud  generally  attained  a  steady  state  condition  as  indicated 
by  the  individual  spectra  measured  over  a  period  of  five  minutes.  Narrower 
droplet  spectra,  predominantly  in  the  radius  range  0.3  to  3.4  pm  were  generated 
by  a  De  Vilbiss  model  65  ultrasonic  nebulizer  with  fixed  frequency  1.35  MHz 
and  a  water  content  output  variable  up  to  6  g  of  water  per  minute.  A  typical 
sequence  of  histogram  giving  the  number  of  drops  per  channel  for  size  range  0 
of  the  CSASP  light  scattering  counter,  is  shown  in  Figure  1.2  for  the  cloud 
nebulizer  generator.  It  can  be  seen  that  the  mode  radius  lies  in  channel  5 
corresponding  to  a  droplet  radius  in  the  range  3.8  to  4.8  pm. 

Water  droplet  clouds  were  also  generated  by  a  comDination  of  a  pair  of  the 
commercial  humidifiers  or  "cool-mist  vapourizers"  at  medium  setting  and  the 
De  Vilbiss  model  65  ultrasonic  nebulizer.  The  cloud  content  output  from  the 
nebulizer  could  be  varied  by  varying  the  amount  of  electrical  power  applied  to 
the  piezoelectric  crystal  in  the  nebulizer.  Averaged  measured  droplet  size 
distributions  using  the  cumbination  of  cloud  generators  described  above  is 
shown  in  Fig.  1.3.  The  cloud  sizes  were  measured  by  the  CSASP  particle  scattering 
counter  and  corrected  in  accordance  with  particle  response  calculations  of 
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Pinnick  and  Auvermann  (1979).  The  number  size  distribution  spectrum  shows 
a  bimodal  structure  and  represents  the  shape  of  the  droplet  size  spectrum 
normally  encountered  at  the  later  stages  of  natural  cloud  development. 

1.4.  Liquid  Water  Content  of  Laboratory  Cloud 

The  prediction  of  atmospheric  extinction  from  a  measureable  fog  or  cloud 
microphysical  parameter  such  as  liquid  water  content  is  of  considerable 
practical  interest.  It  is  important  that  absolute  methods  for  the  measurement 
of  liquid  water  content  are  carefully  evaluated.  Experimental  techniques 
used  to  measure  absolutely  liquid  water  content  of  laboratory  cloud  are 
described  below. 

One  of  the  first  more  reliable  methods  for  the  measurement  of  the  liquid  water 
content  of  fog  was  used  by  Houghton  and  Radford  (1938).  The  amount  of  fog 
water  was  measured  by  passing  a  known  amount  of  fog  laden  air  through  a  series 
of  finely  spaced  thin  wire  screens.  An  impaction  method  is  also  used  here 
whereby  the  cloud  droplets  impact  onto  a  series  of  flannel  filters.  This 
method  was  first  used  by  C.W.  Bruce  (personal  communication)  and  is  referred 
to  by  Bruce  et  al.  (1980).  A  schematic  diagram  of  the  liquid  water  impaction 
device  is  shown  in  Figure  1.4.  The  core  of  the  liquid  water  content  (LWC)  device 
consists  of  an  aluminium  ring  system  which  consists  of  a  series  (usually  4) 
circular  sheets  of  flannel  material  resting  on  a  fine  gauge  metal  perforated 
screen  or  gauze.  A  threaded  collar  is  used  to  firmly  secure  the  filter 
material  when  the  cloudy  air  is  drawn  through  the  device.  An  o-ring  assembly 
is  used  to  ensure  that  no  extraneous  air  is  drawn  into  the  system.  A  compressor 
pump  (P)  together  with  a  calibrated  rotameter  (R)  was  used  to  give  the  volume 
of  cloudy  air  drawn  through  the  filter  assembly  over  a  selected  time  period. 

An  electronic  timer  circuit  in  conjunction  with  a  solenoid  and  a  relay  was 
used  to  actuate  the  aspiration  pumps.  It  also  served  to  remove  a  flap 
positioned  directly  over  the  intake  tube  to  prevent  cloudy  air  from  entering 
the  LWC  device  before  sampling.  The  timer  circuit  was  used  to  preselect  the 
sampling  period  for  the  experimental  measurements.  Subsidiary  comparative 
measurements  were  made  of  the  air  flow  at  the  entrance  and  exit  ports  of  the 
LWC  device  due  to  rela ti vely  smal  1  pressure  drop  in  the  system  and  appropriate 
corrections  were  made  to  the  incoming  airflow  values  as  inferred  from  the 
rotameter  R  readings.  Measurements  using  a  range  of  different  numbers  of 
flannel  filters  indicated  that  4  flannel  filters  were  sufficient  to  capture 
ail  of  the  cloud  water  under  typical  operating  conditions. 
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A  second  identical  LWC  device  was  placed  close  to  the  LWC  impactor  to 
serve  as  a  reference.  Cloudy  air  was  not  aspirated  through  the  reference 
LWC  device,  rather  the  cloud  droplets  were  allowed  to  fall  onto  the  filter 
assembly  over  the  same  sampling  period  as  the  measuring  system.  It  was 
found  that  the  amount  of  cloud  liquid  water  collected  in  the  reference  filter 
was  significant  and  must  be  taken  into  account  in  order  to  take  reliable  and 
accurate  liquid  water  content  measurements.  The  filter  ring  assembly  was 
preweighed  on  an  Oertling  balance  (typical  weight  ~30  g)  and  stored  in  a 
dessicator.  It  was  usual  to  weigh  the  filter,  using  the  Oertling  balance, 
imnediately  after  the  cloud  was  drawn  through  the  LWC  device. 

Another  arrangement  was  also  set  up  to  measure  liquid  water  content  of  the 
cloud  directly.  This  incorporated  a  top  loading  Sartorius  model  1212  MP 
balance.  A  schematic  diagram  of  the  system  is  shown  in  Fig.  1.5.  Cloud  laden 
air  is  drawn  through  a  funnel  shaped  impactor  which  is  packed  with  flannel 
material,  by  means  of  a  compressor  pump.  A  rotameter,  R,  is  also  in  tne  line 
to  measure  the  total  volume  of  cloud  sampled.  A  housing  surrounds  the 
impactor  which  is  mounted  horizontally,  in  order  to  minimise  cloud  from 
entering  the  device  under  conditions  of  no  suction.  The  housing  has  an 
opening  of  diameter  equal  to  the  impactor  tube.  Recordings  of  cloud  water 
content  could  be  made  directly  with  the  top  balance  system.  The  cloud  water 
content  per  unit  volume  was  obtained  from  the  mass  flow  rate  measurements  of  the 
device . 

A  comparison  was  made  between  the  direct  filtration  methods  using  the  vertical 
tube  arrangement  with  its  reference  tube  and  the  top  loading  balance  arrangement 
with  its  horizontal  filter  assembly.  Both  systems  were  adjacent  to  one  another 
in  a  1  m^  chamber  lined  with  water  absorbing  black  material  which  was  pre¬ 
wetted  to  yield  a  water  saturated  environment  for  the  cloud  droplets.  The 
cloud  was  produced  by  a  commercially  available  "cool  mist  vapouriser"  or 
humidifier  whose  output  was  controllable  by  a  range  of  baffle  settings. 

The  results  of  the  LWC  comparison  are  given  in  Table  1.1. 

It  can  be  seen  that  in  general,  agreement  between  the  two  filtration  methods 
is  good.  Agreement  is  somewhat  mitigated  as  the  filter  material  in  the  top- 
balance  assembly  becomes  more  wetted,  resulting  in  an  apparent  loss  of 
collected  cloud  water.  It  is  also  clear  that  the  use  of  a  reference  im¬ 
paction  assembly  is  imperative  when  using  the  direct  filtration  technique 
for  the  measurement  of  liquid  water  content. 
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.5.  Data  Acquisition 

Recording  of  the  radiation  data  was  facilitated  by  means  of  a  digital  panel 
printer  (Datel  Model  DPP-Q7).  The  Data  Acquisition  system  (PDS-200)  was  used 
to  record  and  display  size  distribution  data  in  histogram  form,  as  measured 
by  the  Knollenberg  particle  scattering  counter  -  model  CSASP-100.  A 
Hewlett-Packard  Model  5055A  digital  printer  was  used  to  record  the  particle 
size  channel  data  which  was  displayed  on  the  PDS-200  CRT  screen. 

The  PDP8-E  minicomputer  system  was  linked  successfully  to  the  data  acquisition 
system  (PDS-200)  of  the  particle  scattering  counter.  The  analogue  data  from 
the  radiation  probe  detectors  and  the  electronic  top  balance  was  also  fed  into 
the  PDP8-E  minicomputer  system.  Access  to  the  main  frame  computer  (DEC-20) 
which  contains  the  main  Mie  Lorenz  scattering  computer  routines  (to  be 
described  in  Section  2.2)  was  facilitated  by  a  DEC  SPKL8-JA  terminal  interface 
for  a  9600  baud  rate.  Because  of  the  occurrence  of  a  number  of  failures  in  the 
PDP8-E  mini  computer  (  for  example,  Memory  Bank  connector  failure),  and  the 
serious  disruption  in  the  research  programme  that  would  result  from  a  major 
breakdown  in  the  PDP8-E,  and  the  increasing  difficulty  of  obtaining  component 
parts  it  was  decided  to  opt  for  a  microcomputer  based  data  acquisition  system. 
The  microcomputer  Data  Acquisition  System  is  based  on  a  BBC  microcomputer  (Model 
B).  Two  ports  an  analogue  voltage  imput  port  and  an  8-bit  wide  digital  input/ 
output  port  are  standard  features  of  the  microcomputer.  A  Shugart  SA  300  3£" 
disc  drive  and  an  Epson  F.X-80  line  printer  are  connected  to  the  microcomputer 
providing  permanent  storage  and  hard-copy  display  of  data.  The  monitor  is  a 
"CUB"  Microvitec  2  colour  monitor  which  is  a  medium  resolution  colour  monitor 
and  a  high  resolution  monitor  in  monochrome.  The  microcomputer  can  act  as  a 
computer  terminal  which  is  connected  to  a  DEC  system-20,  a  mainframe  computer 
system.  A  special  Read  Only  Memory  (ROM)  chip  called  a  terminal  emulator  has 
been  installed.  The  program  written  into  the  ROM  was  written  by  the  Computer 
Science  Department  of  Sussex  University,  England.  The  microcomputer  can  emulate, 
or  act  like,  several  types  of  computer  terminals  such  as  the  Digital  Corporation 
VT52.  The  Central  Processing  Unit  or  CPU  is  based  on  a  6502  microprocessor 
and  operates  at  a  clock  frequency  of  2  MHz. 

The  particle  data  Acquisition  system  (PDS-200)  associated  with  the  particle 
counter  was  interfaced  to  the  BBC  microcomputer.  Data  was  transmitted  from 
the  BBC  microcomputer  to  the  DEC  system  20  mainframe  computer,  and  was  facilit¬ 
ated  by  means  of  a  softwate  program  X  READ2  on  the  mainframe  computer.  A 
system  program  called  MLAB  is  used  to  analyse  the  data  which  is  stored  in  the 
user  Directory  on  disk. 
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The  measurements  were  made  in  the  environmental  cloud  housing  of  volume 
1  m  with  the  optical  arrangement  as  shown  in  Fig.  1.1.  Details  of  the 
apparatus  and  ancillary  equipment  used  have  been  already  described  in  the 
previous  sections.  The  vertically  polarised  CC^  laser  beam  traversed  the 
chamber  pathlength  and  emerged  through  a  small  aperture  in  the  exit  wall  of 
the  chamber.  The  aperture  i.i  the  detector  probe  was  reduced  down  to  3  mm 
diameter  in  order  to  minimize  the  entry  of  forward  scattered  radiation  into 
the  detector  (Deepak  and  Box,  1978  a,b). 

A  homogeneous  path  is  usually  assumed  when  the  extinction  coefficient  ae  is 
derived  from  the  Beers -Lambe rt  law  of 

I/I0  =  exp  (-j6L)  (1.8) 

where  I  is  the  incident  radiation  intensity  and  I  is  the  intensity  after 
traversal  of  path  length  L  through  a  (cloud)  medium.  Departures  from  cloud 
homogeneity  will  underestimate  the  inferred  extinction  coefficient. 

A  useful  technique  for  determining  the  extent  of  cloud  homogeneity  has  been 
devised  in  this  laboratory.  The  apparatus  essentially  consists  of  a  rad¬ 
iation  detector  which  is  mounted  inside  the  cloud  housing,  with  the  detector 
being  translated  through  the  cloud  on  a  threaded  rod  assembly.  A  He-Ne 
laser  is  directed  onto  the  detector  which  is  a  UDT  FIL-100V  silicon  photodiode 
operating  in  the  photovoltaic  mode.  The  detector  mounting  is  driven  by  a 
motor  control  module  (whilst  the  scanning  rate  of  the  detector  is  set  by  a 
controlled  oscillator  frequency,  with  5  V  amplitude).  The  direction  of 
rotation  of  the  threaded  rod  can  be  reversed  by  switching  a  100  pF  capacitor 
between  the  motor  inputs.  The  motor  and  control  electronics  are  fan  cooled. 
Scanning  rates  of  between  30  seconds  and  3  minutes  are  readily  achieved  whilst 
mechanical  gears  are  needed  to  extend  this  range. A  narrow  jet  of  cloud  free 
air  is  continuously  directed  across  the  face  of  the  detector  during  a  scan  in 
order  to  prevent  cloud  deposition  on  the  detector  itself. 

A  non-cloud  scan  ensured  a  uniform  response  of  the  detector  along  the  trans¬ 
mission  pa tiu  A  good  representation  of  cloud  homogeneity  along  the  trans¬ 
mission  path  is  shown  in  Fig.  1.6  for  the  cloud  generators  located  in  their 
normal  symmetrical  position  in  the  laboratory  chamber,  where  a  45  degree  slope 
entails  100  per  cent  homogeneity.  Fig.  1.6  shows  the  degree  of  inhomogeneity 
in  the  cloud  over  a  transmitted  path  of  0.7  m  in  terms  of  optical  depth, 

plotted  agaii  path  distance.  Most  cloud  scans  with  this  technique 
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have  yielded  result:,  similar  to  those  in  Fig.  1.6  indicating  a  high  degree  of 
homogeneity  of  cloud  in  the  chamber.  This  is  an  important  finding  for  our 
particular  cloud  chamber  in  that  it  obviates  the  need  for  several  liquid 
water  measurement  devices  to  be  placed  along  the  transmission  path.  It  also 
implies  that  in  general  extinction  measurements  in  the  chamber  do  not  require 
correction  for  inhomogeneity  along  the  transmission  path.  It  should  be 
remarked  that  checks  for  inhomogeneity  are  desirable  particularly  in  chambers 
of  large  dimension.. 

Both  broad  cloud  droplet  size  distribution  using  "cool-riist"  vapourizers  and 
narrow  droplet  size  distributions  using  a  De  Vilbiss  model  65  ultrasonic 
nebulizer  were  used  in  a  1  m  laboratory  chamber.  Simultaneous  measurement 
of  extinction  coefficient  a  ,  m  ^  at  the  CCL  laser  wavelength  A  =  10.591  pm 

e  _3  ^ 

and  liqjid  water  content  W,  gm  using  the  direct  filtration  methods  out¬ 
lined  iri  section  1.4  were  made  using  two  distinct  experimental  procedures: 

The  first  procedure  involved  point  measurements  of  ae  and  W  after  the  cloud 
had  been  allowed  to  reach  a  steady  state  condition.  Point  measurements 

of  aQ  and  W  are  shown  in  Fig.  1.7.  The  error  bars  of  the  liquid  water 
content  parameter  indicates  the  spread  in  liquid  water  using  the  two  filtrat¬ 
ion  methods  described  in  1.4.  Isolated  points  without  error  bars  indicate 
measurements  with  the  vertical  tube  filtration  method  only. 

The  measured  points  reveal  the  ratio  op/W  being  less  than  the  predicted  value 
of  147,  shown  by  the  solid  line  in  Figure  1.7  as  inferred  from  Fq.  (6)  using 
the  analysis  of  Chylek  (1978)  and  Pinnick  et  al  (1979).  The  largest 
deviations  from  prediction  occur  for  the  broader  size  distributions  as  produced 
by  the  cool  mist  vapourizers.  The  value  of  extinction  coefficient  inferred 
from  the  relation 


oe  =  3?c.W/2.\f 

is  overestimated  for  droplet  radius  r  >  14  pm.  Impactor  size  measurements , 
which  will  be  described  in  Section  3  of  this  report  indicate  that  between 
30-35-  of  the  droplet  number  concentration  possess  radii  in  excess  of  14  pm, 
when  produced  by  the  cool  mist  vapourizer  generators.  Accordingly  ce/W  is 
overpredicted  which  is  in  accordance  with  the  measurements  in  Fig.  1.7  which 
are  between  10-20  lower  than  predicted  from  Eq.  (1.9). 

The  second  experimental  procedure  allowed  the  cloud  to  reach  steady  state 
conditions  and  then  the  cloud  generators  were  switched  off  in  order  to  make 
continuous  simultaneous  measurements  of  op  and  W  during  cloud  decay. 

Narrower  size  distributions  produced  from  the  ultrasonic  nebulizer  were  mainly 
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used  for  the  decay  measurements.  Typical  experimental  results  are  shown 

in  Fig.  1.8  and  Fig.  1.9  for  a  liquid  water  content  range  up  to  7  gm  and 

“3 

for  a  lower  liquid  water  content  range  from  2.5gm  down  to  low  values  of  LWC. 

Both  figures  show  good  agreement  with  the  prediction  of  Eq.  (1.9)  shown  by  the 

solid  line  and  they  represent  the  first  set  of  continuous  simultaneous  measure- 

ments  of  a  and  W. 
e 

Thus  the  linear  relationship  as  predicted  by  Chylek  (1978)  between  liquid 
water  content  and  IR  extinction  has  been  experimentally  verified  down  to 
relatively  low  values  of  cloud  water  content  representative  of  natural  cloud 
values.  It  is  also  seen  that  underestimation  of  the  extinction  coefficient 
occurs  for  size  distributions  containing  a  sizeable  proportion  (~30-35/o)  of 
drop  sizes  exceeding  about  14  pm  radius. 
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2 .  EXTINCTION  AND_ BAC KSl A TJX R  IN  WATER  CLOUDS  AT  VISIBLE  WAVELENGTHS 

2.1.  Introduction 

Section  2  describes  numerical  calculations  of  extinction  and  liquid  water 
content  of  cloud  and  fog  at  infrared  and  visible  wavelengths,  based  on  156 
measured  cloud  droplet  size  distributions  and  on  measured  fog  size  distribut¬ 
ions.  Calculations  are  also  presented  of  extinction  and  backscatter  co¬ 
efficient  at  visible  and  near  IR  wavelengths.  The  first  definitive  simultan¬ 
eous  measurements  of  both  oackscatter  coefficient  and  extinction  coefficient 
at  visible  (He-Ne  Laser)  wavelengths  are  described.  Since  a  more  complete 
account  of  this  section  is  contained  within  Appendix  1  (Pinnick,  Jennings 
et  al,  1983)  in  some  detail  and  to  a  lesser  extent  in  Appendix  2  (Jennings, 
1983)  the  more  salient  features  will  only  be  presented  in  this  section. 

2.2.  Mie-Lorenz  Scattering  Computer  Programmes 

Mie-Lorenz  scattering  computer  codes  were  incorporated  successfully  onto  the 
University's  main  frame  DEC-20  computer.  A  DEC  VT  50  VDU  terminal  with 
direct  access  to  the  main  frame  computer  was  acquired.  In  addition  the  Micro- 
vitec  colour  monitor  lined  to  the  BBC  microcomputer  also  served  as  a  terminal 
to  the  DEC-20 


The  Mie  scattering  algorithm  of  Dave  (1968)  used  in  the  Mie-Lorenz  scattering 
code,  which  uses  downward  recurrence  in  the  computation  of  the  complex  function 
A^  is  numerically  stable  but  time  consuming.  The  criterion  for  deciding  between 
upward  recurrence  (is  faster)  and  downward  recurrence  is  presented  by  Wiscombe 
(1980)  and  is  incorporated  into  our  Mie  scattering  algorithms.  When  downward 
recurrence  has  to  be  used,  the  recommendation  of  Wiscombe  (1980)  was  adopted: 
i.e.  downward  recurrence  was  initialized  using  the  Lentz  (1976)  method  rather 
than  Daves  (1969)  method,  requiring  significantly  fewer  number  of  iterations 
to  calculate  An  (mx). 

The  computation  of  scattering  amplitudes  Sj  and  is  achieved  using  Wiscombe 's 
(1980)  algorithm,  and  is  faster  than  the  standard  method  for  particle  response 
and  scattering  cross-section  calculations.  The  small  particle  approximation 
of  Wisccmbe  (1983)  wa  also  inserted  into  the  Mie  scattering  algorithms. 

Values  of  efnciency  factors  using  these  modified  routines  give  agreement  to 
at  least  b  significant  figures  with  values  using  the  original  Dav£  routines. 

We  now  have  facility  to  calculate  reliably  and  accurately: 

(i)  Volume  extinction,  absorption,  scattering  and  backscatter  coefficients 


-  12  - 


for  log-normal,  garona  and  power  law  distributions. 

(i  i  )Efficiency  factors  (in  tabular  and  plot  form)  for  extinction,  absorption, 
scattering  and  backscatter. 

(iii)Forward  scattering  corrections  for  monodispersions  of  particles  for 

extinction  measurements,  in  terms  of  the  path  averaged  correction  factor 
R  (Deepak  and  Box,  1978)  pertinent  to  our  experimental  open  detector 
arrangment. 

(iv) Particle  response  scattering  cross-sections  for  particle  light  scattering 
instruments,  such  as  the  Knollenberg  CSASP-100  particle  counter.  An  example 
of  particle  counter  response  for  the  CSASP-100  instrument  for  water  droplets 
and  NaCl  particles  is  shown  in  Figure  2.1. 

(v)  Volume  extinction,  absorption,  scattering  and  backscatter  coefficients  from 
input  of  data  in  histogram  format.  This  is  of  particular  relevance  to 
computations  using  data  from  the  Knollenberg  particle  scattering  counter. 

2.3.  Middle  IR  Extinction  and  Liquid  Water  Content  in  Cloud 

It  has  been  shown  theoretically  and  verified  numerically  (Chylek,  1978; 

Pinnick  et  al,  1979)  that  an  approximate  linear  relationship  of  the  form  of 
equation  (1.6)  exists  between  middle  infrared  extinction  (around  A  =  11pm) 
and  liquid  water  content  of  fogs.  The  oQ  -  W  relation  has  been  verified 
experimentally  for  laboratory  generated  cloud  (this  work,  Section  1),  Bruce 
et  al  (1980)  and  Gcrtler  and  Steele  (1980).  Natural  cloud  droplets  of  course 
can  be  much  larger  than  those  in  laboratory  cloud  or  fog.  For  this  reason 
we  might  not  expect  relation  (1.6)  to  be  applicable  to  all  clouds,  particularly 
if  droplets  with  radius  r  >  14  pm  dominate  either  extinction  or  liquid  water 
content.  To  investigate  quantitatively  the  magnitude  of  the  error  involved  in 
the  application  of  (1.6)  to  clouds  we  again  made  Mie  calculations  of  the 
extinction  coefficient  and  the  liquid  water  content  for  156  cloud  droplet  size 
distributions  summarized  in  Table  2,  Appendix  1.  The  results  of  these 
calculations  are  compared  to  the  size-  distribution-independent  prediction 
(1.6)  in  Fig.  2.2).  (The  effect  of  gaseous  absorption  is  small  and  has  been 
neglected).  Except  for  cumulonimbus,  nimbostratus,  cumulus  congestus,  orographic 
and  some  stratus  type  clouds  (which  contain  significant  numbers  of  large  (r  > 

14  pm)  droplets)  relation  (1.6)  is  within  a  factor  two  of  the  numerical  results. 
This  comparison  thus  reaffirms  the  conclusion  of  Chylek  (1978)  that  at  between 
10  andll  pm  there  exists  a  nearly-unique  relation  between  extinction  coefficient 
and  liquid  water  content  of  the  form  of  Eq.  (1.6)  for  nonprecipitating  clouds. 
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2.4.  V i s i ble  Ext i  jh; td  on  and  Liqu id  Water  Cent ent  in  Cloud 

A  high  degree  of  correlation  between  visible  extinction  and  liquid  water 
content  of  cloud  or  fog  cannot  be  expected  to  hold  in  general  as  borne  out 
by  the  work  of  Ch#lek  (1978)  and  Pinnick  et  al  (1979).  Chylek  (1978)  has 
shown  that  a  linear  relationship  between  visible  extinction  and  liquid  water 
content,  independent  of  size  distribution  will  only  hold  providing  the  max¬ 
imum  allowable  droplet  radius  does  not  exceed  0.55  pm  -  clearly  at  variance 
for  typical  cloud  and  fog  droplet  sizes.  Indeed  visible  extinction  can  vary 
by  approximately  an  order  of  magnitude  for  a  particular  LWC  value  as  shown 
by  Pinnick  et  al  (1979)  and  reproduced  in  Fig.  2.3. 

Although  Pinnick  et  al  (1978)  found  maximum  radius  values  for  haze  particles 
of  about  2  pm,  some  haze  distribution  will  have  a  significant  proportion 
of  the  size  distribution  less  than  0.55  pm.  In  these  circumstances  linearity 
between  visible  extinction  and  LWC  will  be  expected.  At  near  infrared  wave- 
lrngths  where  the  maximum  allowable  radius  increases  to  1.2  pm  under  which 
conditions  extinction  is  linearly  related  to  LWC,  calculations  by  Pinnick  et 
al  (1978)  indeed  show  linearity  between  and  W  at  wavelength  A  =  1.2pm  for 
the  haze  size  distributions. 

It  was  pointed  out  recently  by  Jennings  (1983)  that  it  is  extremely  unlikely 
that  extinction  will  be  linearly  related  to  number  concentration  for  fog 
(or  cloud)  size  distributions  at  visible  wavelengths. 

2.5.  Extinction  and  Backs catter  in  Water  Clouds  at  Visible  Wavelengths 

Twomey  and  Howell  (1965)  found  that  there  exists  a  relation  (though  not  a 
unique  one)  of  the  form  =  constant  between  backscatter  and 
extinction  coefficient  j  using  Gaussian  and  Poisson  size  distributions. 

Curcio  and  Knes trick  (1958)  made  simultaneous  measuremtns  of  backscatter  and 
extinction  in  cioud  and  found  a  proportional i ty  between  extinction  and  back¬ 
scatter  coefficients  of  the  form  =  for  weathei  conditions  including 

fog,  fog  and  drizzle  and  clear  weather.  However,  there  is  considerable 
leeway  in  determining  the  exponent  in  this  proportional i ty  from  their 
measured  data  (their  ?igun  4).  In  addition  the  effects  of  fog  inhomogeneities 
and  multiple  scatter  contributions  to  the  backscatter  and  transmission 
signals  probably  caused  uncertainty  in  their  data.  Carrier  et  a*i  (1967) 
calculated  extinction  and  backscatter  coefficients  for  eight  model  cloud 
drop  spectra  at  visible  and  middle  IR  wavelengths.  At  0.694  pm  wavelength, 


they  obtain  an  average  value  for  all  8  cloud  models  of  1/22.6  for  the  back- 
scatter  to  extinction  ratio.  More  recently,  Derr  (1980  )  predicts  a  value  of 
1/16.76  for  the  backscatter  to  extinction  ratio  at  i  =  0.694  pm  for  a  set 
of  Deirmendjian  cloud  size  distributions  appropriate  to  cumulus,  continental 
cumulus  and  cumulonimbus  clouds. 

The  backscatter  coefficient  of  a  polydispersion  of  spherical  cloud  drop¬ 
lets  described  by  a  size  distribution  nfr)  and  refractive  index  m  is  given  by 

°b  =  IP  G(m.x)  n(0  dr  (2-1 ) 

where  G(m,x)  is  the  backscatter  gain  defined  as  4ti  times  the  ratio  of  the 
backscatter  differential  cross  section  to  the  geometric  area. 

Cloud  is  dominated  by  droplet  sizes  with  radii  2  pm  <  r<  85  pm  corresponding 
to  size  parameter  12  <  x<  500  at  a  wavelength  A  =  1.06pm.  If  we  can  assume 
that  atmospheric  cloud  droplets  generally  have  slowly  varying  concentrations 
over  radius  intervals  of  the  order  Arssl.6  pm  at  near  IR  wavelengths  (1.06p) 
then  the  Mie  efficiency  factor  for  extinction  Qe  can  be  approximated  (see 
Appendix  1  for  further  details)  by 

Qe  =  2  (1  +  x'2/3)  (2.2) 

The  backscatter  gain  G(m,x)  can  also  be  approximated  by 

G  =  9 { X )  (1  +  Sx2)  (2.3) 

where  g  (A)  is  a  slowly  varying  function  of  wavelength  and  6«1. 


The  use  of  approximations  (2.2)  and  (2.3)  in  equations  (1.1)  and  (2.1) 
leads  to  the  cloud  extinction  coefficient  being  related  to  the  back¬ 
scatter  coefficient  (Pinnick  et  a  1  ,  1983  shown  in  Appendix  1)  and  has  the  form 

-  «•  <  r4/3><  rA  >)  ♦  O (6 2 )1  (2 

<rW>  J 

where  k  :  the  wavenumber  and  <  rn  >  is  the  nth  moment  of  the  droplet  size 
distribution.  If  one  takes  account  of  the  higher  order  terms  in  the  above 
equation,  using  cloud  gamma  size  distribution  models  of  Diermendjian  (1969), 
then  these  higher  order  terms  are  typica  ’y  1Q£  of  the  leading  term  for  non- 
precipi  ta  ti  ng  clouds.  If  errors  of  this  ,  rler  are  acceptable  then  Eq.  (2.4) 
reduces  to  the  simple  form 

Ce  =  T 

where  the  extinction  is  a  linear  function  of  backscatter,  independent  of  drop 
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To  test  the  range  of  validity  of  the  linear  exti nction-backscatter  relation 
(2.5)  for'  cloud  we  c£ilculated  using  Mie  theory  and  indexes  of  refraction  of 
water  given  by  Hale  and  Querry  (1973)  the  extinction  coefficient  according  to 
Eq.(l.l)  and  the  backscatter  coefficient  according  to  Eq.  (2.1)  for  156 
cloud  droplet  size  distributions  measured  in  the  major  cloud  types.  The 
sources  of  these  measurements  together  with  the  range  "of  droplet  sizes  measured 
and  other  pertinent  information  is  listed  in  Table  2,  Appendix  1.  The  main 
sampling  technique  employed  to  obtain  the  cloud  droplet  size  distributions  was 
that  of  impaction  of  droplets  onto  coated  slides  or  replicators  whose  collection 
efficiencies  were  known.  The  practical  lower  limit  for  detection  of  cloud 
droplets  by  the  impaction  technique  is  around  1.5  pm  radius.  The  sole  cloud 
size  determination  by  a  light  scattering  counter  (Ryan,  et  al ,  1972)  was 
calibrated  by  means  of  uniformly  seized  water  droplets.  Only  non-precipitating 
clouds  were  used  in  the  analysis  and  measurements  which  showed  evidence  of 
glaciation  were  excluded. 

The  numerical  calculations  of  extinction  and  backscatter  are  shown  compared 
to  our  linear  relation  (2.5)  in  Fig.  2.4  for  A  =  l.Oym.  (As  an  aside  we  note 
that  Derr's  (1980)  relation  between  extinction  and  backscatter  oe=  18.0  abis 
negligibly  different  from  our  relation  (2.5)  at  this  wavelength.}  For  all 
considered  cloud  size  distributions  the  relation  (2.5)  is  within  50%  of  the 
numerical  results.  Thus  the  numerical  results  suggest  that  cloud  extinction 
coefficients  can  be  inferred  from  measurement  of  the  backscatter  coefficients 
directly  from  Eq.  (2.5),  without  need  to  know  details  of  the  cloud  droplet 
size  distribution.  If  knowledge  of  the  droplet  size  distribution  is  available, 
then  extinction  coefficients  could  be  determined  more  accurately  by  employing 
the  better  approximation  (2.4). 

The  simple  linear  exti nction-backscatter  relation  (2.5)  should  be  particularly 
useful  for  lidar  probing  of  cloud  edges,  where  entrainment  causes  intense 
evaporation  and  severe  distortion  of  the  droplet  spectra,  and  where  as  a  result 
the  spectra  may  not  be  representative  of  the  entire  cloud.  The  associated  wide 
spatial  and  temporal  variability  in  the  droplet  spectra  will  not  prohibit  the 
use  of  our  relation  (2.5)  to  infer  extinction  coefficients  from  lidar  back¬ 
scatter  coefficients,  since  the  relation  is  size-distribution-independent. 

Of  course  cloud  backscatter  coefficients  can  be  determined  from  lidar  return 
signals  (in  a  strai ght-forward  way)  only  in  the  absence  or  neglect  of  multiple 
scattering  contributions  to  the  lidar  signal.  It  follows  that  application  of 
(2.5)  to  obtain  cloud  extinction  coefficients  from  lidar  returns  might  be 
restricted  to  the  edges  of  clouds  where  the  contribution  of  multiply-scattered 
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photons  is  small. 

2.6.  Experimental  Verification  of  the  Extinction-Backscatter  Relation  at  Visible 
Wavelengths. 

The  apparatus  used  tG  measure  the  extinction  and  backscatter  coefficients 
is  shown  schematically  in  Fig.  (2.5).  In  an  attempt  to  determine  experimentally 
the  accuracy  of  the  exti nction-backscatter  relation  (2.5)  we  decided  to  restrict 
our  study  to  laboratory  simulations  of  cloud  (because  of  the  difficulty  in 
making  atmospheric  measurements).  We  generated  in  a  cubical  1.0  m  chamber 
polydispersions  of  cloud  droplets  using  "cool -mist"  vaporizers  and  a  De 
Vilbiss  model  65  ultrasonic  nebulizer  described  in  section  1.3.  Chamber 
saturation  and  cloud  stability  was  maintained  by  soaking  a  black  matt  cloth 
lining  the  chamber  floor  with  water.  The  matt  cloth  also  served  to  reduce 
stray  light  levels.  Droplet  distributions  obtained  with  these  generation 
schemes  were  generally  unimodal  with  a  range  of  drop  sizes  characteristic 
of  clouds.  Interestingly,  the  nebulizer  distributions  were  narrower  than 
those  for  the  cool-mist  vaporizers,  but  they  resulted  in  higher  backscatter 
and  extinction  levels  because  of  their  much  higher  number  concentrations. 

Backscatter  and  extinction  measurements  were  made  simultaneously  on  the 
laboratory  simulated  clouds  employing  a  5.4  itw  He-Ne  laser  source,  a  synchron¬ 
ous  radiometer  (Laser  Precision  Corporation  model  Rk  5100)  for  monitoring 
laser  power,  a  pyroelectric  detector  (LPC  model  RkP-545),  for  measuring  trans¬ 
mitted  laser  power,  and  a  silicon  photodetector  (United  Detector  Technology 
model  FIL-100V)  for  measuring  backscattered  light.  Measurement  of  the  back- 
scattered  light  was  facilitated  with  a  highly  reflective  (>99.3 %)  circular 
mirror  judiciously  positioned  on  the  axis  of  the  laser  beam  so  that  a  small 
(3.16  mm)  hole  drilled  through  its  center  would  admit  the  laser  beam  into  the 
fog  chamber,  and  at  the  same  time  intercept  light  scattered  in  the  near¬ 
backward  direction.  The  mirror  was  tilted  about  9°  (from  the  laser  beam 
direction)  to  reflect  the  backscattered  radiation  onto  the  photodetector.  The 
silicon  photodiode  has  a  responsivity  of  0.40  A  Watt'**  at  X  =  0.6328  pm. 

The  backscatter  signal  was  then  amplified  and  monitored  by  an  accurate  digital 
voltmeter  measureable  to  microvolts.  Test  measurements  of  l^cer  radiation 
were  in  agreement  to  within  1.  using  the  photodiode  detector  and  the  pyro¬ 
electric  detector  probe.  Precautions  were  taken  to  ensure  that  detector 
apera tures  were  small  enough  that  forward  scattering  corrections  to  both 
transmission  and  backscatter  signals  (Deepak  and  Box,  1978)  could  be  neglected. 
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The  extinction  coefficient  oe  (km“^ )  is  determined  from  the  Lambert- Bouguer 
Law  (Eq .1,8)  which  assumes  homogeneity  along  the  cloud  path.  IQ  is  the  measured 
radiation  signal  in  the  absence  of  cloud,  I  is  the  measured  signal  in  the 
presence  of  cloud.  The  backscatter  coefficient  (knf^  sr"^)  is  derived 
from  the  expression 
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where  1^  is  the  measured  backscatter  signal  (in  Watts). 

Ic:  absolute  value  of  the  He-Ne  laser  radiation  (in  Watts). 

5  -1-1 
a^:  the  backscatter  coefficient  (km  sr  ) 

A  :  photodetector  area 

b  :  path  length  between  the  photodetector  and  the  chamber  entrance, 
equal  to  AB+BC  (see  figure  2.5) 

The  effective  photodetector  area  was  reduced  by  8 %  due  to  the  3.16  mm 
diameter  hole  in  the  reflecting  mirror. 


Simultaneous  measurements  of  extinction  and  backscatter  coefficients  were 
recorded  every  ten  seconds  for  each  experimental  run.  The  half  angle  6  sub¬ 
tended  by  the  open  pyroelectric  detector  is  not  more  than  0.1°  and  so  forward 
scattering  corrections  to  extinction  measurements  for  the  droplet  media  are 
less  than  a  few  percent  in  accordance  with  the  calculations  of  Deepak  and  Box 
(1978  a,b).  Measurements  of  radiation  using  a  range  of  entrance  apertures  to 
the  photodiode  detector  indicated  that  multiply  scattered  radiation  did  not 
contribute  to  the  extinction  signal.  The  silicon  photodetector  output  was 
offset  to  a  minimum  signal  level  in  the  absence  of  cloud.  Signal  to  noise 
ratios  >  200  prevailed  for  the  backscatter  measurement. 

Simultaneous  measured  values  of  extinction  and  backscatter  coefficient  made 
at  10  second  intervals  are  presented  in  Figure  2.6,  for  different  experimental 
runs.  The  majority  of  measurements  are  for  the  broader  size  distributions 
produced  by  the  pair  of  humidifiers,  the  narrow  size  distributions  generally 
produced  the  largest  extinction  and  backscatter  values  (data  indicated  by 
O  in  Figure  2.6).  The  ratio  of  extinction  to  backscatter  coefficients  at  the 
He-Ne  wavelength  A  =  0.6328  ym  is  within  20*  of  the  predicted  relation  of 
8 

°e  =  gTX)  where  the  equation  is  a  linear  function  of  backscatter,  independ¬ 

ent  oT  droplet  size,  for  clouds  with  liquid  water  contents  ranging  from  0.05 
up  to  1  gm’^.  Agreement  between  experiment  and  theory  is  best  for  poly¬ 
dispersions  of  droplets  having  a  broad  range  of  sizes  (these  correspond  to  data 
points  in  Fig.  2.6  where  o<  30km  ^sr  ^).The  reason  for  this  can  be  understood 
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by  examining  the  more  exact  backscatter-extinction  relation  (2.4).  For  broad 
distributions  the  higher  order  terms  in  (2.4)  tend  to  cancel  since  they  are  of 
opposite  sign  and  of  comparable  magnitudes,  making  the  zero  order  approximation 
(2.5)  a  good  one.  For  clouds  (generated  with  the  nebulizer  technique)  having 
a  narrow  distribution  of  smaller  droplets  (corresponding  to  data  points  where 
>  30  km  ^sr  ^)  the  agreement  between  experiment  and  theory  is  not  quite  as 
good.  For  these  distributions  the  second  term  in  the  extinction-backscatter 
relation  (2.4)  is  no  longer  concelled  by  the  (negative)  contribution  of  the  term 
of  order  6, rendering  the  linear  (zero  order)  relation  (2.5)  less  accurate. 

The  above  measurements  were  made  solely  for  the  clouds  in  a  steady-state  condition. 
Measurements  were  also  recorded  here  for 

(a)  during  the  initial  build  up  of  cloud, 

(b)  steady-state  cloud  conditions 

(c)  the  initial  decay  interval  of  the  cloud.  The  cloud  size  distributions 

are  generated  by  a  combination  of  two  humidifiers  at  medium  setting  and  the  nebulizer 
whose  size  distribution  is  of  the  form  shown  in  Fig.  1.3.  Simultaneous  measurements 
of  extinction  and  backscatter  coefficients,  recorded  at  10  second  intervals  are 
presented  in  Figures  2.7  and  2.8  for  both  the  cloud  buildup  stage  indicated  by 
"0"  and  the  cloud  decay  stage  "I"  when  the  cloud  generators  were  switched  off. 

In  general,  the  measurements  agree  fairly  closely  (within  about  15%)  with  the 
predicted  size  distribution  independent  relation  -  Eq.  (2.5). 

These  experimental  measurements  -  the  lirst  such  work  that  is  definitive  - 
verifies  the  linearity  of  the  relation 


between  extinction  coefficient  oe  and  backscatter  extinction  at  visible 
wavelengths  (A  =  0.6328  pm)  for  laboratory  generated  fog/cloud  droplet  dis¬ 
tributions.  This  result  suggests  that  visible  or  near-i nfrared  extinction 
coefficients  in  cloud  of  unknown  type  could  be  inferred  from  lidar  backscatter 
measurements  alone,  without  knowledge  of  the  cloud  droplet  size  spectra, 
excluding  complications  thay  may  arise  from  multiple  scattering  contributions  to 
the  lidar  return  signal. 
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3<  EXTINCTION  AND  BACKSCATTER  OF  WATER  CLOUDS  AT  C02  LASER  WAVELENGTHS. 

3.1.  Introduction 

Although  (X^laser  technology,  which  incorporates  CO2  lidar  work, is  now  developing 
rapidly,  efforts  to  extract  attenuation  from  a  CO^  lidar  return  have  to  date  been 
largely  statistical  in  nature.  A  knowledge  of  the  relationship  between  backscatter 
and  attenuation  would  allow  the  determination  of  attenuation  and  backscatter  from 
the  return  signal  of  a  lidar  system  using  the  analytical  inversion  solution  of 
Klett  (1981).  However,  very  little  work  has  been  done  to  date  on  relating  back¬ 
scatter  to  extinction  at  CO^  laser  wavelengths.  Some  of  the  principal  recommend¬ 
ations  which  resulted  from  the  workshop  on  global  large  aerosols,  compiled  by 
Freeman  F.  Hall,  (1983)  Chief  of  the  Doppler  Lidar  program  at  NOAA,  included  the 
necessity  of  CCo  backscatter  and  extinction  measurements  -  both  requiring  experi¬ 
mental  and  theoretical  work. 

We  should  not  necessarily  expect  the  extinction-backscatter  relation  (2.5)  for 
cloud  to  be  applicable  at  all  wavelengths,  since  the  backscatter  gain  cannot 
generally  be  well  approximated  by  slowly  varying  functions  of  size  parameter  at 
all  wavelengths.  To  prove  this  conjecture,  the  extinction  and  backscatter  co¬ 
efficients  for  the  previously  mentioned  156  cloud  size  distributions  were  calcul¬ 
ated  at  several  laser  wavelengths.  An  example  of  the  results  at  the  CO^  laser 
wavelength  A  =  10.6  pm  (Fig.  7  Appendix  1)  show  that  for  a  particular  backscatter 
coefficient  the  extinction  varies  by  an  order  of  magnitude  for  different  size 
distributions  of  droplets. 

Despite  the  scatter  in  the  calculated  points  of  o^  of  Fig.  7  (Appendix  1),  a 
closer  examination  of  the  results  reveal  that  the  broader  drop  size  distributions 
within  the  156  cloud  size  distributions  investigated  shows  that  the  extinction  to 
backscatter  ratio  attains  relatively  constant  values. 

The  normalized  backscatter  cross-section  (efficiency  factor)  Q  (m  ,x)  for  water 
is  plotted  in  Fig.  3.1  as  a  function  of  size  parameter  x  using  an  updated  Mie- 
Lorenz  scattering  computer  code  using  Wiscombe's  (1980)  algorithm,  details  of 
which  are  described  more  fully  in  Section  2.2.  It  is  seen  that  the  normalized 
backscatter  cross-section  oscillates  about  a  constant  value  for  x  >  6  (radius 
r  >  10  pm  at  A  =  10.591  pm).  In  view  of  the  fact  that  calculations  show  that 
the  extinction  cross-section  is  constant  for  increasing  x  for  x  >  8  at  A  = 

10.591  pm,  we  can  predict  that  exti ncti on/backscatter  will  De  largely  independent 
of  x  also  for  x  >  6-3. 


a 
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Calculations  indeed  show  that  the  extinction  to  backscatter  ratio  ae/ab  at 
A  -  10.591  ym  (P20  C02  laser  wavelength  line)  oscillates  about  a  constant 
value  for  water  droplet  radius  >  8.5  ym  as  shown  in  Fig.  3.2.  We  predict 
that  broad  cloud  drop  size  distributions  will  give  rise  to  relatively  constant 
values  of  extinction  to  backscatter  ratios  at  CO2  laser  wavelengths. 

For  size  parameter  x<  8  the  efficiency  factor  for  extinction  Qe  can  be 
apDroximated  by  the  linear  relation  (1.4),  Qe  (m,x)  =  c(A).  x  at  C02  laser 
wavelengths.  The  backscatter  gain  G(m,x)  is  a  rapidly  varying  function  of  x 
(see  Fig.  3.1)  over  the  same  size  parameter  range.  It  follows  that  for  x <  8 
(droplet  radius  <  14  ym)  at  C02  laser  wavelengths,  the  extinction  to  backscatter 
ratio  is  a  sensitive  function  of  drop  size  as  is  borne  out  by  the  numerical 
calculations  of  extinction  and  backscatter  for  natural  cloud,  shown  in  Fiq.  7 
(Appendix  1). 
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3.2.  Experimental  Measurement  of  Backscatter  Coefficient  and  Extinction  Coefficient 
at  CC>2  Laser  Wavelengths. 

( a )  Experimental  Apparatus  and  Techniques 

Backscatter  and  extinction  measurements  were  made  on  laboratory  simulated 
clouds  employing  the  model  941S  CO2  laser  source,  a  synchronous  radiometer 
(model  Rk51 00)  for  monitoring  laser  power,  a  pyroelectric  detector  (LPC  model 
RkP-545)  for  measuring  the  transmitted  laser  power  and  an  Infrared  Associates 
Inc.  Mercury  Cadmium  Telluride  HgCdTe  liquid  nitrogen  cooled  8-13pm  radiation 
photoconducti ve  detector  (Model  No.  HCT-90)  for  measuring  the  backscattered 
light.  Attempt  were  made  to  use  a  similar  technique  described  in  Section 
2.6  for  visible  wavelengths  whereby  a  central  hole  drilled  through  a  metalized 
mirror  surface  would  admit  the  laser  beam  through  the  cloud  chamber  and  at 
the  same  time  intercept  the  backscattered  radiation  onto  the  HgCdTe  detector. 
However  it  was  not  possible  to  use  this  technique  because  of  the  occurrence 
of  directly  reflected  radiation  from  the  internal  walls  of  the  hole  onto  the 
sensitive  detector.  The  detector  was  placed  close  to  the  CC^  laser  axis  in 
the  backward  direction,  at  a  distance  not  more  than  5.9  cm  from  the  axial 
direction,  thus  subtending  an  angle  not  less  than  177.4°  with  the  forward 
direction  at  the  chamber  centre.  The  HgCdTe  detector  was  mounted  on  an  x-y 
translation  stage  with  vertical  height  adjustment. 

The  evaluated  detector  possessed  a  metal  dewar  assembly  which  contained  the 
necessary  liquid  nitrogen  for  cooling  the  detector.  The  HCT-90  detector  had 
a  sensing  area  of  1  mm  squared.  The  signal  from  the  detector  was  fed  into  a 
PPA-15-1  pre-amplifier,  then  into  the  input  of  a  EG&G  Brookdeal  Model  9501 
lock-in  amplifier.  A  EG&G  Brookdeal  model  9479  light  chopper,  operating 
normally  at  1000  Hz  served  as  reference  to  the  lock-in  amplifier.  The  amplifier 
output  from  the  lock-in  amplifier  was  recorded  on  a  Philips  model  PM8100  chart 
recorder.  Precautions  were  taken  to  reduce  stray  reflected  IR  laser  radiation 
from  entering  the  detector  during  experimental  runs.  The  detector  aperture 
on  the  pyroelectric  probe  RkP-545  was  maintained  at  3  mm  to  minimise  forward 
scattering  corrections  to  the  transmission  signal  (Deepak  and  Box,  1978). 

The  water  clouds  were  produced  using  "cool -mist"  vapourisers  on  their  own 
and  in  combination  with  a  De  Vilbiss  model  65  ultrasonic  nebulizer  (see 
section  1.3).  Typical  droplet  size  spectra  obtained  using  a  two-stage 
impactor  (Garland,  1971)  mounted  in  a  wind  tunnel  are  shown  in  Fig.  3.3, 
using  the  "cool-mist"  humidifiers  at  medium  (upper  figure)  and  high  settings. 
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particle  scattering  counter)  to  be  measured.  The  droplets  which  impacted 
onto  gelatine  slides  were  viewed  and  measured  microscopically  with  the  aid 
of  a  Porton  graticule. 

( b )  Measurement  of  Backsca tter  Coefficient  and  Extincti on  Coefficient  at 
CO9  Laser  Wavelengths. 

Due  to  the  necessarily  high  sensitivity  of  the  HgCdTe  detector,  it  was  not 
possible  to  make  simultaneous  measurement  of  extinction  and  backscatter  in 
the  chamber  cloud.  (This  is  because  the  backscattered  or  reflected  signal 
from  the  walls  of  the  pyroelectric  detector  probe  at  the  exit  end  of  the 
chamber  greatly  exceeded  the  true  signal).  The  experimental  procedure  in¬ 
volved  firstly  the  measurement  of  the  initial  CC^  laser  intensity  (IQ) 
before  entering  the  chamber  followed  by  measurement  of  the  backscatter 
signal , interspersed  with  extinction  measurements.  Both  backscatter  and 
extinction  measurements  were  made  for  steady-state  cloud  conditions  and 
during  cloud  decay.  The  backscatter  coefficient  (km”^  Sr~^)  is  derived 
from  equation  (2.6).  The  absolute  value  of  the  CO2  laser  radiation,  in 
Watts,  was  determined  from  the  pyroelectric  RkP-545  probe  which  was  calib¬ 
rated  by  means  of  an  Electro  Optical  Industries  Inc.  Blackbody  Source, 
model  WS  153,  of  one  inch  conical  cavity  with  temperature  controller  prov¬ 
iding  a  temperature  range  from  50  to  1000°C. 

The  measured  backscatter  signal  with  output  from  the  lock-in  amplifier  was 
calibrated  by  directing  a  known  low  CO^  laser  radiation  signal  onto  the 
HgCdTe  detector.  The  low  level  CO2  laser  radiation  signals  were  monitored 
on  the  model  Rk  5100  radiometer  operating  in  synchronization  with  a  CTX-32 
external  chopper, 

Siqnal  to  noise  ratios  of  at  least  >  10  (and  generally  much  in  excess  of  10) 
prevailed  for  the  backscatter  measurements.  Measurements  of  backscatter 
coefficient  and  extinction  coefficient  at  a  =  10.591  pm  (unless  otherwise  stated) 
are  presented  in  Fiq.  3.4.  Cool-mist  vapourizers  were  used  here  to  produce  the 
laboratory  cloud.  Measured  values  taken  in  steady  st'te  cloud  conditions,  are 
indicated  by  the  points  "8"  in  Fig.  3.4.  The  experimental  points  "I"  represent 
measurements  of  backscatter  and  extinction  coefficients  made  during  the  growth 
staqe  of  the  cloud  up  to  the  steady  state  value. 

Most  measured  values  of  the  extinction  to  beckscatter  ratio  are  above 

the  lower  bound  line  representing  a  value  of  350  in  Fig.  3.4.  This  lower 
bound  value  also  represents  the  minimum  values  obtained  for  oe/o^  using  numerical 


calculations  of  backscatter  and  extinction  coefficient  based  on  156  measu 
droplet  size  distributions  (see  Fig.  7,  Appendix  1).  The  vast  majority  of 

measurements  for  the  cool-mist  humidifier  cloud  yield  values  of  extinction  to 
backscatter  ratio  from  350  up  to  550,  as  shown  in  Fig.  3.4.  These  ratio  values 
reflect  the  fact  that  the  cloud  droplet  size  parameters  were  predominantly  <  than 
about  6  (radius <  about  10  micrometres)  ,  a  region  where  the  backscatter  cross- 
section  possesses  much  greater  values  than  average  (Fig.  3.1). 


Fig.  (3.2)  predicts  that  the  extinction  to  backscatter  ratio  is  a  sensitive 
function  of  drop  size  for  droplet  radius<  10  pm.  This  prediction  is  borne 
out  by  the  upper  measurements  in  Fig.  3.4  which  were  performed  on  a  different 
type  of  cloud  :  produced  by  a  combination  of  the  ultrasonic  nebulizer  with  two 
humidifiers.  The  narrower  size  distributions  produced  the  larges  backscatter 
and  extinction  signals  because  of  their  much  higher  concentrations.  It  is  clear 
that  the  measured  extinction  to  backscatter  ratio  is  sensitive  to  the  cloud  size 
distribution  over  the  range  of  droplet  size  used. 

Measurements  were  generally  made  at  the  CO^  laser  P20  line  (X  =  10.591  pm). 
However  some  measurements  were  made  at  the  R  18  line  (A  =  10.261  pm)  for  cloud 
produced  by  the  combination  of  nebulizer  and  two  humidifiers.  These  are  shown 
in  Fiq.  (3.4)  and  are  indicated  by  "X".  They  yield  larger  c^/a^  ratio  values, 
principally  because  the  extinction  per  unit  mass  is  greater  at  A  =  10.261  pm  than 
at  A=  10.591  pm  (see  Pinnick  et  al ,  1979  for  details). 


Conclusions 

The  measurements  described  in  this  section  represent  the  first  definitive 
measurements  of  backscatter  and  extinction  coefficient  in  laboratory  cloud  at 
CO^  laser  wavelengths.  The  measurements  are  in  reasonably  good  agreement  with 
numerical  predictions  on  cloud  possessing  narrow  type  size  distributions  as  are 
used  here. 

The  measurements  yield  relatively  constant  values  for  extinction  to  backscatter 
ratio  for  the  same  common  form  of  size  distribution.  The  ratio  values  vary 
between  350  and  550  in  these  experiments.  It  is  also  clear  that  the  measurements 
show  sensitivity  of  0^/0^  to  size  distribution  when  a  different  shaped  size 
distribution  is  used. 
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Fig.  1.3  Measured  laboratory-generated  cloud  droplet 

size  distributions  used  to  examine  the  accuracy 
of  the  extinction-backscatter  relation.  The 
droplet  generators  consisted  of  two  humidifiers 
(at  medium  setting)  and  a  De  Vilbiss  Model  65 


cloud  at  wavelength  X  -  10.591  v  • 
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1.8  Continuous  simultaneous  measurements  of  ae  and  W  for  laboratory  cloud 
at  X  =  10.591  pm.  The  lower  solid  line  is  a  best  fit  to  the  data. 

The  upper  straight  line  is  that  predicted  by  the  Ch^lek  relation. 
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Fig.  2.1  Mie  theory  response  calculations  for  the  Knollenberg  CSASP  particl 
counter  for  NaCl  particles  (lighter  curve)  and  water  droplets. 


SCATTERING  CROSS  SECTION  BETWEEN  4.00  DEGRES  AND  22.00  DEGRES 
N  =  I .332000,  K  =  -1 .470000E-08  ,  VAVELEN6TH  =  6.328080E-0I  MICRONS 


EXTINCTION  COEFFICIENT 


Fig.  2.2  Volume  extinction  coefficient  at  A=1 0. 6  ym  versus  liquid  water 
content  for  156  cloud  droplet  size  distribution  measurements  of 
cumulus  and  stratus  clouds. 
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Fig.  2.4  Volume  extinction  coefficient  versus  volume  backscatter  coefficient 
at  wavelength  X  =  1.06pm  for  156  droplet  size  distributions 
measured  for  the  major  cloud  types. 
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Schematic  diagram  of  the  experimental  system  used  to  measure  extinction  and  backscatter 


Fig.  2.o  Comparison  of  measured  cloud  backscatter  and  extinction 
coefficients  with  the  theoretical  relation  (2.5). 
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Fig.  3.3  Histogram  size  distribution  measurements  of 
laboratory  cloud  using  a  two  stage  impactor. 
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Backscatter  and  Extinction  in  Water  Clouds 

r.  G.  Pinnick,1  S.  G.  Jennings,2  Petr  Chylek,3  Chris  Ham,4  and  W.  T.  Grandy,  Jr.3 

An  approximate  relation  between  the  volume  extinction  coefficient  at  and  backscatter  coefficient  eb  of 
atmospheric  cloud  at  visible  and  near-infrared  wavelengths  is  derived.  The  relation  is  only  weakly 
dependent  on  the  size  distribution  of  droplets  and  has  the  form  om,'o%  «  (8n/g){  1  +  k~  ,/J«/'4,,)/<r:» 

—  5Ck1«/,4>/</,*>)  4-  k*'J«r4'J)<r*>/<r2><r*»]}  where  the  extinction  efficiency  is  approximated  by  a 
complcx-angular-momentum-theory  result  and  the  parameters  g  and  <5  are  determined  by  approximating 
a  running  mean  of  the  backscatter  gain  by  G(x)  -  j(l  +  <5x:)  (x  is  droplet  size  parameter  and  <5  «  l).  k  is 
the  wave  number,  and  <r">  is  the  nth  moment  of  the  droplet  size  distribution.  To  zero  order  the  relation 
is  linear  and  independent  of  the  droplet  size  distribution  a,  [8rr/g<A)3<y*  where  ^(i)  is  a  slowly  varying 
function  of  wavelength.  At  a  wavelength  X  »  1.06  >tm  the  relation  is  <r,  (km'1) «»  18.2  ab  (km-1  sr-1). 
Predictions  made  with  this  simple  zero-order  approximation  are  in  good  agreement  (within  50%)  with 
Mic  calculations  of  extinction  and  backscatter  coefficients  based  on  J56  measurements  of  cloud  droplet 
jpecira  in  cumulus  and  stratus  type  clouds.  The  linear  tr#  —  ob  relation  is  also  in  agreement  with 
extinction  and  backscatter  measurements  made  on  laboratory- generated  fog  droplet  distributions.  The 
relation  suggests  that  visible  or  near-infrared  extinction  coefficients  in  cloud  of  unknown  type  could  be 
inferred  from  lidar  backscatter  measurements  alone,  without  knowledge  of  the  doud  droplet  size  spectra, 
barring  complications  arising  from  multiple  scattering  contributions  to  the  lidar  return. 


1.  Introduction 

A  possible  relation  between  atmospheric  backscatter  and 
extinction  (or  visibility)  has  been  of  interest  to  scientists  for 
more  than  2  decades.  Curcio  and  Knesrriek  [1958]  found  ex¬ 
perimentally  by  using  a  white  light  source  a  dehnite  cone- 
la  tion  between  backscatter  and  transmission  measurements 
for  atmospheric  conditions  ranging  from  relatively  dear  wea¬ 
ther  to  fog  to  drizzle.  Twomey  and  Howell  [1965]  were  able  to 
at  least  partially  explain  Currio's  findings  by  showing  that  for 
polydispersions  of  water  droplets  the  ratio  of  calculated  parti¬ 
cle  backscatter  to  extinction  for  a  white  light  source  is  oi)ly  a 
weak  function  of  droplet  size  distribution.  They  used  Gaus¬ 
sian.  Poisson,  reversed  Poisson,  and  bimodal  distributions 
characteristic  of  atmospheric  fog  and  haze  in  their  numerical 
study.  Additional  experimental  evidence  for  correlation  of 
backscatter  and  extinction  for  white  light  sources  was  provid¬ 
ed  by  Vogi  [1968], 

With  the  coming  of  the  laser  used  in  the  lidar  technique  and 
its  application  to  the  remote  measurement  of  cloud  [Co//is, 
1965;  Schotland  et  ai,  1971;  Plass  and  Kaitawar,  1971;  Zuev 
and  Bolin,  1972;  Platt ,  1973;  Derr  et  al .,  1976;  and  others] 
there  was  additional  interest  in  the  question  of  a  possible 
relation  between  backscatter  and  extinction,  but  in  this  case 
for  monochromatic  sources.  Although  Silverman  and  Sprague 
[1970]  did  not  intend  to  utilize  a  laser  source  in  their  instru¬ 
ment  for  single-ended  measurements  of  visibility,  their  feasibil¬ 
ity  studies  are  for  monochromatic  sources.  They  calculated  the 
backscatter  coefficient  ab  and  the  extinction  coefficient  cr,  for 
polydispersions  of  water  droplets  having  gamma-type  size  dis¬ 
tributions  considered  to  be  characteristic  of  common  fog  and 
cloud  types.  Their  numerical  studies  indicate  an  approximate 
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relation  between  backscatter  and  extinction  of  the  form  a ,  ~ 
trt0  90.  Similar  numerical  investigations,  but  specifically  for  the 
lidar  application,  have  been  made  by  Derr  [1980].  Den- 
showed  the  backscatter  to  extinction  ratio  for  ruby  (X  *»  0.694 
//m)  and  Nd-YAG  ( X  1.06  /zm)  sources  is  a  slowly  varying 
function  of  drop  size.  Like  Silverman,  Derr  calculated  back¬ 
scatter  to  extinction  ratios  averaged  over  various  gamma  size 
distributions  representative  of  cumulus  and  cumulonimbus 
cloud  and  from  these  results  estimated  approximate  linear 
relations  between  backscatter  and  extinction  coefficients. 

Up  to  this  point,  then,  investigators  have  found  either  ex¬ 
perimentally  or  through  numerical  studies  that  there  should 
exist  at  least  an  approximate  relation  between  backscatter  and 
extinction  coefficients  in  water  clouds,  but  the  underlying 
reason  for  the  relation  is  unclear. 

In  this  paper  we  re-examine  from  a  theoretical  standpoint, 
in  light  of  recent  complex-angular-momentum  theory  predic¬ 
tions  [Nussenzveig,  1969;  Khare  and  Nussenzveig,  \911a\Nuss- 
enzveig  and  Wiscombe,  1980],  the  question  of  a  possible  rela¬ 
tion  between  backscatter  and  extinction,  and  the  universality 
such  a  relation  might  have. 

We  derive,  using  a  complex-angular-momentum  theory  ap¬ 
proximation  for  the  extinction  efficiency,  and  a  simple  curve- 
fit  approximation  for  the  running  mean  of  the  backscatter 
gain,  a  relation  between  extinction  and  backscatter  that  to 
zero  order  is  independent  of  droplet  size  distribution;  higher- 
order  terms  in  our  solution  take  the  size  distribution  depen¬ 
dence  into  account.  These  terms  are  worked  out  explicitly  for 
gamma-type  size  distributions  characteristic  of  cloud  and  are 
found  to  contribute  on  the  order  of  10%  of  the  leading  term. 
Because  of  the  simple  form  of  the  solution  which  neglects 
higher-order  terms  (wherein  extinction  is  related  linearly  to 
backscatter  with  no  "size  distribution  dependence),  we  con¬ 
centrate  on  investigating  the  universality  of  this  (zero  order) 
solution.  We  compare  the  zero  order  solution  (which  at 
X  -  1.06  pm  is  a,  (km-1)  »  18.2  ab  (km-1  sr-1))  to  both  Mie 
calculations  (for  156  measured  droplet  distributions)  and 
measurements  of  backscatter  and  extinction  made  on  labora¬ 
tory  generated  clouds.  Agreement  is  in  both  cases  generally 
within  50%.  Our  experimental  verification  of  this  linear  a, 
—  relation  is  the  first  such  work  that  is  definitive. 
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Also  shown  are  results  for  the  lowest  order  CAM  theory  terms  used 
m  the  derivation  of  the  backscatter-extmction  relation  (17). 


relation  exists  at  an  infrared  laser  (10.6  pm)  wavelength,  where 
the  backscatter  to  extinction  ratio  can  vary  by  more  than  an 
order  of  magnitude  with  the  form  of  the  cloud  droplet  size 
distribution.  We  also  show  that  cloud  liquid  water  content  for 
clouds  of  unknown  droplet  si2e  distribution  cannot  be  inferred 
from  visible,  infrared,  or  near-millimeter  backscatter  measure¬ 
ments  alone.  Finally,  we  demonstrate  that  previously  derived 
relations  between  infrared  (a  =  10.6  pm)  extinction  and  liquid 
water  content  in  fog.  and  between  infrared  (a  =  3-8  pm)  ab¬ 
sorption  and  liquid  water  content  in  fog  [CAy/efc.  1978; 
Pinmck  er  al. ,  1979]  can  also  be  applied  to  most  clouds. 

2.  Backscatter  and  Extinction  in  Cloud 

The  volume  extinction  and  backscatter  coefficients  a,  and 
e,  of  a  polydispersion  of  spherical  cloud  droplets  char¬ 
acterized  by  a  size  distribution  n(r )  and  refractive  index  m  are 
given  by 

a,  « 

J  sr2Gfl(r)  dr  (2) 

where  Q,[m,  ,t)  is  the  Mic  efficiency  factor  for  extinction  for  a 
particle  with  refractive  index  m  and  size  parameter  x  *=  Inr/X, 
and  Gim,  x)  is  the  backscatter  gain  defined  as  4x  rimes  the 
ratio  of  the  backscatter  differential  cross  section  to  the  geo¬ 
metric  area.  Because  of  the  complex  size  dependence  of  Q4  and 
G,  if  wc  are  to  find  an  approximate  relation  between  extinc¬ 
tion  and  backscatter  coefficients,  we  will  no  doubt  have  to 
resort '  to  some  approximate  expressions  for  the  Mie  ef¬ 
ficiencies  Q '  and  G  in  equations  (1)  and  (2). 

3.  Approximation  for  the  Extinction  Efficiency 
Extinction  in  cloud  is  dominated  by  droplets  with  radii  2 
|jm  <  r  <  85  p m,  corresponding  to  size  parameters 
12  <  x  <  500  at  a  wavelength  X  »  1.06  n m.  These  relatively 
large  size  parameter  values  suggest  asymptotic  expansions  for 
the  Mic  extinction  and  backscatter  efficiencies,  and  just  such  a 
theory  has  been  developed  by  Nusscnzveig  [1969],  who  ana¬ 
lytically  continued  the  Mie  series  into  the  complex  angular 
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momentum  plane.  The  resulting  analytic  expressions  for  the 
scattering  functions,  for  large  x.  are  given  in  terms  of  the  poles 
and  saddle  points  contributing  to  the  scattering  amplitudes  in 
the  complex  plane  of  summation  index.  The  compfox-angular- 
momentum  (CAM)  theory  expression  for  the  extinction  ef¬ 
ficiency  has  been  worked  out  by  Nussenzveig  and  Wiscombe 
[1980]: 

n  2 

e.  -  2  t  1.992.x-:'3  -  — ; - sin  [2<m  -  1)JC]X-  1 

(m  —  IKm  i) 

-0.715x“*'3  +  0(x~5'3)  (3) 

A  plot  of  this  approximate  expression  for  Q4  (Figure  1)  shows 
the  familiar  damped  oscillating  behavior  with  period  Ax  =  */ 
(m  —  1)  s  9.5.  If  we  can  assume  atmospheric  cloud  droplets 
generally  have  slowly  varying  concentrations  over  correspond¬ 
ing  radius  intervals  (Ar  =  1.6  /im  at  X  *=  1.06  /im),  then,  since 
we  are  concerned  only  with  integrals  of  Q„  over  cloud  size 
distributions,  wc  can  neglect  the  oscillatory  sin  term  in  (3).  We 
can  also  neglect  the  terms  0(.x~40)  and  Oix'5'3)  since  for 
.x  >  10  they  are  less  than  ~  10%  of  the  second  term  0(x"2'3). 
Given  these  constraints  we  can  approximate  Q ,  by 

<?.  =  2(1  +  (4) 

as  shown  in  Figure  1. 


4.  Approximation  for  the  Backscatter  Gain 

The  backscatter  cross  section  (and  backscatter  gain)  can 
also  be  approximated  by  using  complcx-angular-momentum 
theory.  Within  the  framework  of  CAM  theory  the  single¬ 
particle  backscatter  cross  section 
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(5) 


is  evaluated  by  making  a  Debye  expansion  of  the  scattering 
matrix  [iVussenrueig,  1979] 


SM'  *)  -  Z  S)JX’  *)  J  -  1.  2  (6) 

r-o 

and  identifying  dominant  terms.  In  this  expansion  the  pth 
term  is  associated  with  [p  —  1)  internal  reflections  at  the  dro¬ 
plet  surface,  except  for  p  «  0,  which  is  associated  with  direct 
reflection.  According  to  the  computer  evaluations  of  Khare 
and  Nuzzenzveig  [Nussenzveig,  1979]  for  water  droplets  with 
index  of  refraction  given  by 

n  -  [cos  (1  la/48)]- 1  -  1.33007  (7) 

the  following  four- term  approximation  to  the  amplitudes  ac¬ 
counts  for  80-90%  of  the  backscattered  intensity 

^(x,  K)  =  SjJ""  +  Sjj ,< '»  +  +  Shs*  (8) 

These  arise  from  the  (p  =*  2)  residue  series,  the  lOth-order  rain¬ 
bow,  and  the  geometrical  optics  rays  al  p  «=  0,  p  ®  2.  The 
intensity  computed  from  (8),  ISjj1,  gives  rise  to  interference 
among  the  four  terms,  resulting  in  quasiperiodic  oscillations 
with  periods  Ax  *=  0.41,  0.83,  1.1,  14  [Nussenzveig,  1979]. 
These  periods  are  evident  in  Mie  calculations  of  the  back¬ 
scatter  gain  in  Figure  2. 

Since  our  interest  is  in  backscattering  from  a  polydispersion 
of  many  cloud  droplets  of  different  size,  rather  than  back- 
scattering  from  single  droplets,  we  are  motivated  to  consider 
averages  of  the  intensity.  If  one  averages  over  the  largest 
period.  Ax  *=  14,  the  intensity  in  the  backward  direction  is  well 
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Fig.  2.  Mic  calculations  or  the  back  scatter  gun  Gt.x)  for  water  droplets  at  a  wavelength  X  -  1.06  pm  (refractive  index 
m  -  1.325  -  5  *  10 ■'/)  calculated  for  size  parameters  around  x  =  300  with  resolution  tn  size  parameter  Jkr  **  0.01.  The 
quasipenodic  structure  with  periods  Ax  =  0.41,  0.83,  14  are  apparent. 


approximated  by  the  sum  of  the  squares  of  the  amplitudes  on 
the  right-hand  side  of  (8).  Moreover,  the  geometrical  optics 
terms  are  CXx)  and  make  small  contributions  in  the  range 
20  <  x  <  500.  Thus,  for  our  purposes  (water  drops  in  cloud), 

IS,!*.  n)|:  -  |S;,.""’(x.  *>|;  +  |S,.11I'V  x)|»  (9) 

and  only  the  p  «  2  residue  scries  and  p  =  11  {10th-order  rain¬ 
bow;  Debye  terms  of  the  scattering  amplitude  need  be  evalu¬ 
ated. 

The  following  asymptotic  expressions  for  these  terms  can  be 
extracted  from  /Chare's  [1976]  thesis 

Tenth-Order  Rainbow  (/>  =  //) 

5t  n  {R)(x.t7)  =  TT7TTfrv?exp  (kA  ){/>,/■  —  polt)AHsr  iCl) 

-  k"  l!\q0c  -  (10) 

Where  k  =  2.x,  the  parameters  A,  p0r,  p0y ,  q0*t  and  q0 M  are 
determined  by  saddle  point  evaluations  in  the  complex  angu¬ 
lar  momentum  plane  [Khurc  and  Aiusscnveig.  1977h).  and  At  is 
the  Airy  function  [ Abramowit :  and  Stegun .  1966],  Evaluation 
of  these  parameters  leads  to  the  following  expression  for  the 
lOth-ordcr  rainbow  term 

|5»  tt) j*  =  1 7 5.99 a'10'3 [9. 8  x  10~1Mi(-0.0(mix2;;>)|2 

+  1.37  x  10_  5x“2'3Mf( -0.0041  lx2'3)j:]  (11) 

Residue  Series  (p  •=  2) 

Once  again  we  extract  from  Kharc’s  thesis  the  leading  be¬ 
havior  of  this  term: 

S,  x.  n)  =  —  0.24276x4'3cxp  (-0.42074x,/3) 

•exp  [i(2mx-  -  rr/6  -f  x(:)]  (12) 

where  Cj  is  a  small  purely  real  angle.  Finally,  squaring  this 
amplitude  gives 

\St'2irn\x.  zr)l2  -  0.05S93xB''3exp(-0.84I5x,,3)  (13) 


The  CAM  theory  result  for  the  back  scatter  gain,  calculated 
according  to  equations  (5),  (9),  (II),  and  (13) 

G{x)  =  704x4'3[9.8  x  10' 7Mr(- 0.0041  lx2'3)  |2 

+  1.37  x  10' 5x“2/s|i4r(— 0.0041  lx2,3)|2] 

+  0.236x2'3exp  (-0.842x,/3)  (14) 

are  shown  compared  to  the  exact  Mie  result  (averaged  over 
Ax  «=  0.83)  in  Figure  3. 

Although  the  CAM  theory  prediction  has  the  correct  form 
(it  is  slowly  varying  over  the-range  20  <  x  <  500,  but  never¬ 
theless  steadily  increases  for  x  >  200),  it  is  too  low  by  an 
additive  constant  of  about  1.2.  (As  an  aside  we  note  that  the 
geometrical  optics  ray  contribution  to  the  backscatter  gain 
has  the  correct  form  to  account  for  this  discrepancy,  but  it  is 
still  small  compared  to  the  terms  in  (14).)  We  offer  the  follow¬ 
ing  reasons  for  the  underprediction  by  CAM  theory:  (1)  We 
have  retained  only  leading  order  corrections  in  the  asymptotic 
expansions  of  the  p  -  2,  11  terms  in  the  Debye  series,  and  (2) 
there  is  no  apparent  fundamental  reason  for  the  existence  of  a 
unique  functional  relation  between  backscattering  and  extinc¬ 
tion  coefficients,  for  they  correspond  to  entirely  different 
physical  processes.  That  the  theory  nevertheless  provides  quft!» 
itativdy  correct  verification  of  our  expectations  is  remarkable, 
and  thus  the  connection  seems  worthy  of  further  analysis.  We 
leave  this  further  analysis  for  future  work. 

Because  of  our  failure  to  reproduce  the  correct  magnitude 
result  with  CAM  theory,  we  instead  propose  a  simple  curve-fit 
approximation  to  the  backscatter  gain.  Since  our  interest  is  in 
lidar  backscattering  from  a  polydispersion  of 'droplets  (rather 
than  from  single  droplets)  we  again  neglect  the  oscillatory 
component  evident  in  the  exact  Mie  results  and  approximate 
the  gain  by 


G(x)  -  g(X)([  +  ax2)  (15) 


where  g(X)  is  a  slowly  varying  function  of  wavelength  and 
<5  «  1.  This  approximate  expression  is  shown  compared  to  the 
exact  Mi:  results  in  Figure  3. 
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TABLE  2.  Cloud  Size  Distribution  Measurements 


Source 

Cloud  Type 

Range  of  Droplet 
Sizes  Measured 
(radii  in  pm) 

Number  of  Drop 
Size  Distribution 
Measurements 

aufm  Kampe  and 

Cumulus  Congestus 

1.5-92 

1 

Weickmann  [1952] 

Battan  and  Rtitan 
[1957] 

Cumulus 

Cumulus  Congestus 

Tropical  Cumulus 

1.75-58 

5 

Diem  [1943] 

Cumulus 

Cumulus  Congestus 
Stratocumulus 

Altostratus 

Nimbostratus 

Stratus 

1-42 

6 

Durbin  [1959] 

Cumulus 

0.75-30 

22 

Eagan  ei  al. 

[1974] 

Stratocumulus 

Cumulus 

1.25-15 

12 

Fiizgerald  [  1972] 

Continental  Cumulus 
Maritime  Cumulus 

3.5-11.5 

7 

Fitzgerald 

and  Spvers-Duran 
[1973] 

Jiusto  [1967] 

Cumulus 

Stratocumulus 

1.75-10.5 

4 

Maritime  Cumulus 
Orographic  Stratocumulus 

1.5-24 

4 

Ryan  et  al.  [1972] 

Continental  Cumulus 
Maritime  Stratus 

Maritime  Cumulus 

2-42 

33 

Singleton  and 

Stratus 

1.5-62 

17 

Smith  [1960] 

Spvers-Duran 

[1972] 

Altostratus 

Altocumulus 

2.5-24 

8 

Squires  [1958] 

Orographic 

Tradewind  Cumulus 
Continental  Cumulus 

15-82 

10 

Warner  [1969.  1973a] 

Cumulus 

1.25-24 

20 

Warner  [19736] 

Maritime  Cumulus 

1.5-13 

4 

Weickmann  and 

aufm  Kampe  [1953] 

Cumulus  Congestus 
Cumulonimbus 

15-100 

3 

capitation  cloud  model  C.6.  Although  the  series  appears  to 
-  diverging  for  the  C.6  model,  this  is  not  the  case.  In  fact,  for 
recipitating  clouds  the  higher-order  terms  decrease  rapidly  in 
magnitude  and  further  partially  cancel  the  term  CM<5)  in  equa- 
ons  (17)  and  (19),  rendering  the  leading  term  more  dominant, 
his  suggests  that  if  ~20%  errors  are  acceptable,  equation 
;7)  has  a  particularly  simple  form 


here  the  extinction  is  a  linear  function  of  backscatter,  inde¬ 
pendent  of  droplet  size. 

Of  course,  the  argument  that  higher-order  terms  in  (19)  are 
mall  assumes  cloud  droplet  distributions  are  broad  (so  that 
he  oscillatory  terms  in  Qt  and  G  can  be  neglected)  and  further 
hat  they  have  a  single-mode  gamma-tvpe  size  distribution. 
This  latter  assumption  in  particular  is  not  very  realistic  for  all 
loud  types.  In  the  next  section  we  investigate  the  range  of 
aliditv  of  the  linear  extinction-backscatter  relation  (20)  for 
nore  realistic  (measured)  cloud  distributions. 

6.  Numerical  Verification  of  the 
Extinction-Backscatter  Relation 
To  test  the  range  of  validity  of  the  linear  extinction- 
backscatter  relation  (20)  for  cloud  we  calculated,  using  Mie 


Fig.  4.  Volume  extinction  coefficient  versus  volume  backscatter 
coefficient  at  a  wavelength  X  *  1.06  pm  for  156  droplet  size  distri¬ 
butions  measured  in  the  major  cloud  types:  Cu,  denotes  cumulus, 
cumulus  congestus,  continental  cumulus,  maritime  cumulus,  tropical 
cumulus,  altocumulus,  and  tradewind  cumulus;  Sc,  stratocumulus; 
Cb,  cumulonimbus;  Or,  orographic;  St,  stratus;  As.  altostratus;  Ns, 
nimbostratus;  and  Mst  denotes  maritime  stratus.  The  results  are  in 
good  agreement  with  the  size-distribution-independent  prediction  (20) 
(shown  by  the  solid  straight  line)  relating  extinction  uniquely  to  back- 
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r  ^  Tvpica!  Uboratorv-eenerated  cloud  droplet  distributions 
.  JdVmve^iMie  the  accuracC  of  the  cxtmction-backscattcr  relation 
POT  wasurements  Were  made  with  a  model  CSASP-100  light 
scattering  counter  (manufactured  by  Particle  Measurement  Systems, 
Inc  Boulder.  Colo.)  calibrated  for  water  droplets  according  to  the 
prescription  of  Pmmck  and  Auvermann  [1979]. 


theory  and  indexes  of  refraction  of  water  given  by  -Hale  and 
Querrv  [1973],  the  extinction  coefficient  according  to  equation 
(I)  and  the  backscattcr  coefficient  according  to  equation  (2)  for 
156  cloud  droplet  size  distributions  measured  in  the  major 
cloud  types.  The  sources  of  these  measurements  together  with 
the  ranee  of  droplet  sizes  measured  and  other  pertinent  infor¬ 
mation  is  listed  in  Table  2.  The  main  sampling  technique  em- 
ploved  to  obtain  rhe  cloud  droplet  size  distributions  was  that 
of  impaction  of  droplets  onto  coated  slides  cr  replicators 
whose  collection  efficiencies  were  known.  The  practical  lower 
limit  for  detection  of  cloud  droplets  by  the  impaction  tech¬ 
nique  ts  around  1.5  ^im  radius.  The  sole  cloud  size  determi¬ 
nation  by  a  light  scattering  counter  [Ryan  et  ai,  1972]  was 
calibrated  by  means  of  uniformly  sized  water  droplets.  Only 
nonprecipitating  clouds  were  used  in  the  analysis  and 
measurements  which  showed  evidence  of  glaciation  were  ex¬ 
cluded. 

The  numerical  calculations  of  extinction  and  backscatter 
are  shown  compared  to  our  linear  relation  (20)  in  Figure  4  for 
/.  «  1.06  fim.  (As  an  aside  we  note  that  Derrs  [1980]  relation 
between  extinction  and  backscatter  cre  -  18.0  crb  is  negligibly 
different  from  our  relation  (20)  at  this  wavelength.)  For  all 
considered  cloud  size  distributions,  the  relation  (20)  is  within 
50";.  of  the  numerical  results.  If  the  measured  droplet  distri¬ 
butions  were  extrapolated  beyond  the  range  of  measured  sizes 
(imposed  by  the  instrumentation),  then  the  calculated  extinc¬ 
tion  and  backscatter  coefficients  would  probably  fall  even 
closer  to  our  relation  (20).  The  reason  is  that  extrapolation  of 
a  distribution  will  allow  more  chance  for  cancellation  of  the 
oscillatory  terms  in  Qt  and  G,  which  are  neglected  in  deri¬ 
vation  of  (20).  Thus,  the  numerical  results  suggest  that  cloud 
extinction  coefficients  can  be  inferred  from  measurement  of 
the  backscatter  coefficients  directly  from  (20)  without  need  to 
know  details  of  the  cloud  droplet  size  distribution.  If  knowl¬ 
edge  of  the  droplet  size  distribution  is  available,  then  extinc¬ 
tion  coefficients  could  be  determined  more  accurately  by  em¬ 
ploying  the  better  approximation  (16). 

The  stmple  linear  extinction-backscatter  relation  (20)  should 
be  particularly  useful  for  lidar  nrobinc  of  cloud  edges,  where 


entrainment  causes  intense  evaporation  and  severe  distortion 
of  the  droplet  spectra  and  where,  as  a  result  the  spectra,  may 
not  be  representative  of  the  entire  cloud.  The  associated  wide 
spatial  and  temporal  variability  in  the  droplet  spectra  will  not 
prohibit  the  use  of  our  relation  (20)  to  infer  extinction  coef¬ 
ficients  from  lidar  backscattcr  coefficients  since  the  relation  is 
size  distribution  independent. 

Of  course,  backscatter  coefficients  within  cloud  can  be  de¬ 
termined  from  lidar  return  signals  only  so  far  as  radiation 
from  the  transmitter  can  penetrate  the  cloud  and  further  be 
backscattered  with  sufficient  intensity  so  as  to  be  detectable  by 
the  receiver.  This  penetration  depth  is  on  the  order  of  several 
tens  of  meters  for  heavy  cloud  having  liquid  water  content  of  a 
few  tenths  of  a  gram  per  cubic  meter.  Probing  the  interior  of 
clouds  has.  in  addition  to  the  limitations  posed  by  loss  of 
signal  due  to  attenuation,  complications  arising  from  multiple 
scattering  contributions  to  the  lidar  return.  Monte  Carlo  cal¬ 
culations  of  Kunkel  and  Weinman  [1976]  and  Platt  [1981] 
show  these  contributions  are  appreciable  for  cloud  optical 
depths  greater  than  about  0.5,  depending  on  the  solid  angle 
subtended  by  the  lidar  receiver.  It  is  beyond  the  scope  of  this 
paper  to  discuss  methods  for  taking  multiple  scattering  contri¬ 
butions  into  account.  We  only  point  out  there  are  compli¬ 
cations  in  the  interpretation  of  lidar  data  that  cannot  be  over¬ 
come  simply  by  using  knowledge  of  our  relation  between 
backscatter  and  extinction420).  On  the  other  hand,  this  rela¬ 
tion  should  simplify  the  interpretation  of  lidar  data,  even 
when  multiple  scattering  contributions  to  the  lidar  return 
signal  are  important. 

7.  Experimental  Verification  of  the 
Extinction-backscatter  Relation 

It  would  be  desirable  to  compare  the  extinction-backscatter 
relation  (20)  to  direct  measurements  of  these  quantities.  The 
only  known  simultaneous  measurements  of  backscatter  and 
extinction  in  atmospheric  cloud  are  by  Curcio  and  Knestrick 
[1958].  They  found  empirically  a  proportionality  between  ex¬ 
tinction  and  backscatter  coefficients  of  the  form  ot  **•  a*0,64 
for  weather  conditions  including  fog,  fog  and  drizzle,  and  clear 
weather.  However,  there  is  considerable  leeway  in  determining 
the  exponent  in  this  proportionality  from  their  measured  data 
(their  Figure  4).  In  addition,  the  effects  of  fog  inhomogeneities 
and  multiple  scattering  contributions  to  both  the  backscatter 
and  the  transmission  signals  are  potential  uncertainties  in 
comparing  our  relation  (20)  with  their  data. 

In  an  attempt  to  determine  experimentally  the  accuracy  of 
the  extinction-backscatter  relation  (20)  we  decided  to  restrict 
our  study  to  laboratory  simulations  of  cloud  (because  of  the 
difficulty  in  making  atmospheric  measurements).  We  gener¬ 
ated  in  a  cubical  1.0  m3  chamber  polydispersions  of  cloud 
droplets  using  ‘cool-mist’  vaporizers  and  a  De  Vilbiss  model 
65  ultrasonic  nebulizer.  Chamber  saturation  and  cloud  stabili¬ 
ty  was  maintained  by  soaking  a  black  matt  cloth  lining  the 
chamber  floor  w’ith  water.  The  matt  cloth  also  served  to 
reduce  stray  light  levels.  Droplet  distributions  obtained  w’ith 
these  generation  schemes  were  generally  unimodal  (Figure  5) 
with  a  range  of  drop  sizes  characteristic  of  clouds.  Interest¬ 
ingly,  the  nebulizer  distributions  were  narrower  than  those  for 
the  ccol-mist  vaporizers,  but  they  resulted  in  higher  back¬ 
scatter  and  extinction  levels  because  of  their  much  higher 
number  concentrations. 

Backscatter  and  extinction  measurements  were  made  simul¬ 
taneously  on  the  laboratory  simulated  clouds  employing  a  5.4 
mw  He-Ne  laser  source,  a  synchronous  radiometer  (Laser  Pre- 
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„  Corporalion  model  Kk  MOW  lot  . . .  laser 

,  »  pvroelectric  dc.coto,  (U’C  ntt.de!  Kkl’-NM,  for 

Laso'r.ns  transmitted  laser  power,  and  a  silicon  photo- 
Z'Z,  (United  Detector  Technology  mode!  FIU-IOOV)  for 
measuring  backscallered  light.  Measurement  of  the  back- 
jMercd  light  was  facilitated  with  a  highly  reflective  (> 

W  5-.)  circular  minor  judiciously  positioned  on  the  axis  of 
the  laser  beam  so  that  a  small  (3.16  mm)  hole  drilled  through 
us  center  would  admit  the  laser  beam  into  the  fog  chamber. 
Ind  at  the  same  time  intercept  light  scattered  in  the  near- 
tuckward  direction.  The  mirror  was  tilted  about  9°  (from  the 
User  beam  direction)  to  reflect  the  backscattcrcd  radiation 
onto  the  photodetcctor.  Precautions  were  taken  to  ensure  that 
detector  aperatures  were  small  enough  that  forward  scattering 
corrections  to  both  transmission  and  backscatter  signals 
[Dttpak  and  Box .  1978]  could  be  neglected. 

The  results  of  these  measurements  are  summarized  in 
Figure  6  where  the  ratio  of  extinction  to  backscatter  coef¬ 
ficients  at  the  He-Ne  wavelength  X  —  0.6338  pm  is  within  20% 
of  that  predicted  by  our  size-distribution-independent  relation 
(20)  for  clouds  with  liquid  water  contents  ranging  from  0.05  to ..., 
10  %  m"J.  (In  this  comparison  we  have  used  ;/(a  *  0.6943 
=•  1.42  rather  than  a  value  of  if  for  X  =  0.6328  pm.  but  the 
ditTcrenct  is  estimated  to  be  smaii.)  Agreement  between  experi¬ 
ment  and  theory  is  best  for  polydispersions  of  droplets  having 
a  broad  range  of  sizes  {these  correspond  to  data  points  in 
Figure  6  where  o„  <  30  km  ~  1  sr  ’).  The  reason  for  this  can  be 
understood  by  examining  the  more  exact  backscatter- 
extinction  relation  (17).  For  broad  distributions  the  higher- 
order  terms  in  (17)  tend  to  cancel  since  they  arc  of  opposite 
sign  and  of  comparable  magnitudes,  making  the  zero-order 
approximation  (20)  a  good  one.  For  clouds  {generated  with 
the  nebulizer  technique)  having  a  narrow  distribution  of 


BACKSCATTER  COEFFICIENT  UrtiV) 


Fig.  7.  Same  as  Figure  4  except  for  the  CO,  laser  wavelength 
a  =  10.6  pm.  At  this  longer  wavelength  the  extinction  is  no  longer 
approximately  related  to  backscatter,  as  the  extinction  coefficient  for 
a  particular  backscatter  coefficient  varies  by  about  an  order  of  mag¬ 
nitude  with  the  droplet  size  distribution. 

tinction  and  backscatter  are  in  overall  better  agreement  with 
theory  (relation  (20))  than  our  numerical  results  in  section  6, 
which  are  based  on  measured  atmospheric  cloud  droplet  dis¬ 
tributions.  Part  of  the  explanation  for  this  finding  has  to  do 
with  the  limited  scope  of  our  laboratory  studies.  While  we 
tried  to  generate  droplet  distributions  characteristic  of  atmo¬ 
spheric  cloud,  the  range  of  size  distributions  was  nowhere  near 
that  for  those  used  in  the  numerical  study.  Had  we  been  able 
to  generate  larger  droplets  and  narrower  distributions  in  the 
laboratory,  no  doubt  the  agreement  between  measurements 
and  the  linear  backscattcr-extinction  relation  (20)  would  not 
have  been  as  good. 


smaller  droplets  (corresponding  to  data  points  where  ob  >  30 
km  ~  1  sr*5)  the  agreement  between  experiment  and  theory  is 
not  quite  as  good.  For  these  distributions  our  neglect  of  the 
oscillatory  contributions  to  the  extinction  and  backscatter  ef¬ 
ficiencies  render  the  linear  (zero  order)  relation  (20)  less  accu¬ 
rate. 

Nonetheless,  it  is  noteworthy  that  our  measurements  of  ex- 
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Fig  6.  Comparison  of  measured  cloud  backscatter  and  extinction 
coefficients  (points)  with  the  theoretical  relation  (20).  The  backscatter 
and  extinction  coefficient  measurements  were  made  on  laboratory- 
generated  clouds  having  a  variety  of  size  distributions  (two  examples 
uf  which  arc  shown  in  Figure  5)  and  liquid  contents  ranging  from  0.05 
to  t.O  g  m-3.  Clouds  corresponding  to  the  square  data  points  were 
generated  wiih  the  ultrasonic  nebulizer  technique;  all  others  were 
generated  with  cool-mist  vaporizers. 


8.  Extinction  and  Backscatter  at  Middle  Infrared 
AND  MtLUMETER  WAVELENGTHS 

We  should  not  necessarily  expect  the  extinctiotvbdcUcafter 
relation  (20)  for  cloud  to  be  applicable  at  all  wavelength*, 
since  the  backscatter  gain  cannot  generally  be  well  approxi¬ 
mated  by  slowly  varying  functions  of  size  parameter  at  all 
wavelengths.  To  prove  this  conjecture,  we  calculated  the  ex¬ 
tinction  and  backscatter  coefficients  for  the  previously  men¬ 
tioned  156  cloud  size  distributions  at  several  laser  wavc- 
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Fig.  8.  Volume  backscatter  coefficient  at  a  wavelength  4  -  1.06 
/jm  versus  liquid  water  content  for  156  measured  droplet  size  distri¬ 
butions  of  cumulus  and  stratus  type  clouds.  The  results  show  doud 
liquid  water  content  is  not  uniquely  related  to  the  backscatter  coef¬ 
ficient  irrespective  of  cloud  type. 
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Fig.  9  Same  as  Figure  8,  except  for  X  *  10.6  pm 

Wirth5  An  example  of  the  results  at  the  CO,  laser  wavelength 
r.  10  6  pm  (Figure  7)  show  that  for  a  particular  backscatter 
Coefficient  the  extinction  varies  by  an  order  of  magnitude  for 
different  sire  distributions  of  droplets.  (Our  neglect  of  gaseous 
absorption,  which  is  at  most  0.4  km-  at  this  wavelength. 
doesTot  significantly  affect  the  results  in  Figure  7.)  Therefore. 
,  CO  lidarmeasurement  could  not  be  used  (by  itself,  without 
constraints  on  spatial  variability  of  drop  size  distnbut.ons)  to 
6cdu«  infrared  (A  -  10.6  pm)  extinction  in  cloud.  (As  an  aside 
we  note  that  if  the  form  of  the  distribution  of  droplets  is 
spatially  invariant  the  method  of  Klett  [1980]  can  be  used  to 
infer  extinction  profiles  from  backscatter  measurements,  re- 
eardless  of  the  particular  form  or  the  distribution.) 

Neither  can  an  approximate  extinction-backscatter  relation 
be  expected  to  hold  at  millimeter  wavelengths  as  can  be  seen 
from  examining  (1)  and  (2),  knowing  that  the  Rayleigh  ap¬ 
proximation  holds  for  which  G(x)  ~  x*  and  Q.(x)  ~  x.  Our 
numerical  results  based  on  the  156  drop  distributions  (not 
shown  here)  bear  out  this  conclusion. 


backscatter  coefficient  only  through  a  parameter  that  depends 
on  droplet  size  distribution.  To  obtain  a  quantitative  measure 
of  this  size  distribution  dependence,  we  again  performed  Mie 
calculations  of  the  backscatter  coefficient  by  using  cquution  (2) 
and  the  liquid  water  content  by  using  equation  (21)  for  the 
previously  mentioned  156  cloud  size  distributions.  The  results 
are  presented  in  Figure  8  and  show  that  for  a  particular  back¬ 
scatter  coefficient  the  cloud  liquid  water  content  can  vary  by 
more  than  a  factor  or  10  with  the  droplet  size  distribution. 

Similar  investigations  of  a  possible  relation  between  cloud 
liquid  water  content  and  backscatter  coefficient  at  other  infra¬ 
red,  visible,  and  near-millimeter  laser  wavelengths  A  -  0.55, 
0.694,  3.8,  10.6,  1364  pm  (220  GHz),  2143  pm  (140  GHz),  and 
3192  pm  (94  GHz)  show  again  that  no  unambiguous  relations 
exist;  and  further  that  for  a  fixed  backscatter  coefficient  at 
these  other  wavelengths  theUloud  liquid  water  content  is  gen- 
erally  an  even  more  sensitive  function  of  the  droplet  size  distri¬ 
bution.  An  example  of  these  results  at  X  -  10.6  pm  is  shown  in 
Figure  9. 

We  can,  therefore,  conclude  that  for  cloud  with  droplets  of 
unknown  size  distribution  a  determination  of  liquid  water 
content  cannot  be  made  solely  from  a  single-wavelength  lidar 
measurement. 

10.  Extinction  and  Liquid  Water  Content  in  Cloud 

It  has  previously  been  shown  theoretically  and  verified  nu¬ 
merically  {ChyUk,  1978;  Pinnick  et  aL,  1979]  that  approxi¬ 
mate  relationships  exist  between  infrared  extinction  (around 
X  m  1 1  /im)  and  liquid  water  content  of  fogs,  and  between 
infrared  absorption  (around  X  3.8  /im,  9.5  pm)  and  liquid 
water  content  of  fogs.  The  relationships  are  linear  of  the  form 
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(23) 

2  Xp 

(24) 
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9.  Backscatter  and  Liquid  Water  Content  in  Cloud 

Having  been  encouraged  by  the  success  of  the  extinction- 
backscatter  relation  (20)  at  visible  and  near-infrared  wave¬ 
lengths.  we  extended  our  investigation  to  see  if  a  similar  rela¬ 
tion  might  exist  between  cloud  liquid  water  content  and  back¬ 
scatter  coefficient;  the  motivation  of  course  being  the  prospect 
of  utilizing  lidar  for  remote  measurement  of  cloud  liquid  water 
content. 

The  liquid  water  content  \V  of  cloud  with  droplet  size  dis¬ 
tribution  n(r)  is  given  by 

W  -  p  Jy  rMr)  Jr  (21) 

where  p  is  the  density  of  water. 

We  already  know  the  backscatter  gain  for  polydispersions 
of  droplets  at  visible  and  near-IR  wavelengths  can  be  approxi¬ 
mated  by  G(x)  *s  g{\  +  i5x2).  Hence,  there  can  be  no  size- 
distribution-independent  relation  between  liquid  water  con¬ 
tent  and  backscatter  coefficient  at  this  wavelength  as  the  ratio 
of  these  quantities  (after  expanding  in  powers  of 

(22, 

3 9  L<r  >  <r  ><>  >  J 

contains  in  the  leading  term  the  ratio  of  the  third  to  second 
moments  of  the  droplet  size  distribution.  In  other  words,  at 
X  «  1.06  pm  the  liquid  water  content  of  cloud  is  related  to  the 


where  ot  and  am  are  the  extinction  and  absorption  coefficients 
at  the  wavelength  A,  W  the  fog  liquid  water  content,  and  the 
parameters  c  and  c'  are  equal  to  the  slopes  of  straight  lines 
that  approximate  the  Mic  extinction  and  absorption  efficiency 
curves  by  Q£x,  2)  =»  c{X)x  and  Q,(x,  A)  -  c\X)x.  The  <7,  -  W 
relation  has  been  verified  experimentally  for  laboratory- 


Fig.  10.  Volume  extinction  coefficient  at  a  wavelength  X  -  10.6 
pm  versus  liquid  water  content  for  156  cloud  droplet  size  distribution 
measurements  of  cumulus  aod  stratus  douds.  Except  for  cumu¬ 
lonimbus,  nimbostratus,  cumulus  congest  us,  orographic,  and  some 
stratus  cloud  types  the  results  are  dose  to  the  equation  (23)  prediction 
(shown  by  the  straight  line)  relating  infrared  extinction  coeflitient  to 
liquid  water  content.  f  r  *  •  .  — -•  • 
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Fig  11.  Volume  absorption  cociTiciem  at  a  wavelength  X  **  3.8 

versus  liquid  waier  content  for  156  cloud  droplet  size  distri¬ 
butions  of  cumulus  and  stratus  clouds  For  most  cloud  types  the 
results  are  close  to  the  equation  04)  prediction  (shown  by  the  straight 
line)  relating  cloud  infrared  absorption  unambiguously  to  cloud 
liquid  water  content. 

generated  fogs  by  Center  and  Steele  [1930]  and  Bruce  et  al. 
[1930],  The  success  of  relations  123)  and  (24)  depend  on  the 
fact  that  fog  droplets  have  radii  predominately  less  than 
r  S  14  um  [ Ckvlek ,  1978,  Pinnick  et  al 1979]. 

Cloud  droplets  of  course  can  be  much  larger  than  those  in 
fog.  For  this  reason  we  might  not  expect  relations  (23)  and  (24) 
to  be  applicable  to  all  clouds,  particularly  if  droplets  with 
radius  r  >  14  ^m  dominate  either  extinction,  absorption,  or 
liquid  water  content  To  investigate  quantitatively  the  mag¬ 
nitude  of  the  error  involved  in  the  application  of  (23)  and  (24) 
to  clouds,  we  again  made  Mie  calculations  of  the  extinction 
and  absorption  coefficients  and  the  liquid  water  content  for 
the  previously  considered  156  cloud  droplet  size  distributions 
summarized  in  Table  2.  The  results  of  these  calculations  are 
compared  to  the  sizc-disiribution-independcni  predictions  (23) 
and  (24)  in  Figures  10  and  11.  (The  effect  of  ga:.eous  absorp¬ 
tion  is  small  and  has  been  neglected.)  Except  for  cumu¬ 
lonimbus.  nimbostratus,  cumulus  congestus.  orographic,  and 
some  stratus  type  clouds  [which  contain  significant  number  of 
large  ( r  >  14  /im)  droplets)  the  relations  (23)  and  (24)  are 
within  a  factor  of  2  of  the  numerical  results.  This  comparison 
thus  reaffirms  the  conclusion  of  Chylek  [1978]  that  at  /.  =  11 
jrm  there  exists  a  nearly  unique  relation  betwen  extinction 
coefficient  and  liquid  water  content  of  the  form  of  (23)  for 
nonprecipitating  clouds.  It  also  suggests  the  absorption-liquid 
content  relation  (24)  can  be  applied  to  most  clouds  without 
regard  to  their  type  or  character  of  their  droplet  size  distri¬ 
bution. 

11.  Conclusion 

For  all  types  of  atmospheric  clouds  consisting  of  spherical 
water  droplets,  an  approximate  relation  between  their  extinc¬ 
tion  and  backscatter  coefficients  at  visible  and  near-infrared 
wavelengths  has  been  derived  To  zero  order  the  relation  is 
independent  of  cloud  droplet  size  distribution.  The  relation 
should  enable  the  determination  of  cloud  extinction  coefficient 
(or  total  droplet  surface  area)  solely  from  lidar  return  signals, 
providing  the  contribution  of  multiply-scattered  photons  to 
the  lidar  return  can  be  neglected.  However,  no  size- 
distribution-indcoendcnt  relation  exists  between  cloud  liquid 
water  content  and  backscatter  coefficient  at  visible,  infrared. 


or  near-millimeter  wavelengths,  suggesting  that  single- 
wavelength  lidar  (by  itself)  cannot  be  used  to  remotely  mea¬ 
sure  cloud  liquid  water  content  for  clouds  of  unknown  type. 
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Recently,  Lenham  and  Clay1  obtained  high  correlation 
for  a  linear  relationship  between  a  measured  extinction 
coefficient  in  fog  at  wavelength  X  =  0.53  /im  and  liquid  water 
content.  The  liquid  water  content  values  were  derived  from 
number-size  distribution  histograms,  which  in  turn  were  in¬ 
ferred  from  transmission  measurements  at  several  wave¬ 
lengths  using  an  inversion  algorithm.2  They  also  found  a  high 
correlation  between  the  inferred  number  concentration  in  the 
largest  size  category  (6.4~12.8-pm  radius)  and  the  extinction 
coefficient. 

The  purpose  of  this  Letter  is  to  point  out  that  high  corre¬ 
lation  between  visible  extinction  and  liquid  water  content  for 
fog  cannot  be  expected  to  hold  in  general.  In  addition,  the 
correlation  of  number  concentration  with  extinction  is  ex¬ 
amined,  and  the  conclusion  is  drawn  that  it  is  very  unlikely 
in  genera!  for  extinction  to  be  linearly  related  to  number 
concentration  for  fog  size  distributions  at  visible  wave¬ 
lengths. 

The  volume  extinction  coefficient  ae  is  given  by 

c,  =  (  atil(m,x)nir)dr,  (1) 

Jo 

where  n(r)  is  the  droplet  size  distribution.  The  single -panicle 
extinction  cross  section  aexl{m,x)  is  a  function  of  the  droplet 
complex  index  of  refraction  m  and  the  size  parameter  x  = 
2irr/\,  the  ratio  of  the  droplet  circumference  to  the  wave¬ 
length  X. 

Chylek3  has  showm  that  the  extinction  ar  is  linearly  related 
to  liquid  water  content  W  by 

(2* 

2pA 

where  p  is  the  density  of  water,  and  the  coefficient  c  is  equal 
to  the  slope  of  a  straight  line  that  approximates  the  normal¬ 
ized  extinction  cross  section  Qext(m,x)  by 

Q,lt(m pt )  -  c[m)x  (3) 

Chylek1  has  shown  that  for  X  =  0.5  pm,  Q^xi  is  proportional 
to  x  with  c  =  0.61  providing  fog  droplets  have  a  size  parameter 
of  <6  5.  A  similar  result  is  obtained  here  for  X  =  0.53  pm 
(used  by  Lenham  and  Clay1).  This  corresponds  to  a  maxi¬ 
mum  allowalm  radius  of  0.55  pm  if  a  linear  relationship  be¬ 
tween  extinction  and  liquid  water  content  iLWO  is  to  hold 
independent  of  size  distribution.  It  is  clear  since  the  fog 
droplet  size  distributions  of  Lenham  and  Clay1  possess  radii 
greatly  in  excess  of  0.55  pm  that  no  unique  relation  should 


Table  I.  Mean  Extinction  Cross  Sections  per  Droplet  lor  the  Listed 
Channel  Widths  at  Wavelength  0.53-pm;  Index  of  Refraction  m:  1.3338- 
1.45  X  10~* 


Mean  extinction  cross  section 

Channel 

Radius  interval 

per  particle 

(pm) 

(cm2) 

1 

0.2-0.4 

7.655  X  10-9 

2 

0.4-0.8 

3.859  X  10-8 

3 

0.8-1.6 

1.131  X  10-7 

4 

1. 6-3.2 

4.180  X  10"7 

5 

3.2-6.4 

1.603  X  10-6 

6 

6.4-12.8 

6.254  X  10-6 

12.8-25.6 

2.464  X  10-5 

exist  between  extinction  and  LWC  at  visible  wavelengths. 
This  result  also  applies  to  fog  size  distributions  in  general. 

A  Mie  scattering  codo  which  incorporates  Wiscombe's4 
recent  algorithms  is  used  to  compute  the  mean  extinction  cross 
sections  per  particle  for  the  histogram  channel  widths  used 
by  Lenham  and  Clay1  and  for  an  additional  channel  width  of 
12.6-25.6  pm.  The  interpolated  value  of  the  complex  re¬ 
fractive  index  of  water5  at  X  =  0.53  pm  was  used.  The  com¬ 
puted  mean  extinction  cros»  sections  are  given  in  Table  I. 
The  contribution  to  the  extinction  coefficient  by  an  individual 
channel  is  obtained  from  the  product  of  mean  extinction  cross 
section  and  droplet  number  concentration  of  that  channel 
(assuming  the  number  distribution  is  uniform  over  the 
channel). 

The  mean  extinction  cross  section  increases  approximately 
as  the  square  of  the  radius  for  a  large  size  parameter,  as  can 
be  seen  in  Table  I  for  the  last  four  radius  intervals.  For  an 
equal  number  concentration  of  cm-3  over  their  radius  inter¬ 
vals.  channel  6  dominates  (76%  contribution)  the  extinction 
from  these  channels.  An  examination  of  the  differentia! 
number  concentration  of  cm-3  pm-1  for  the  three  largest  size 
intervals  (4, 5,  and  6  in  Table  I)  in  Fig.  2  of  Lenham  and  Clay1 
reveals  that  the  total  extinction  is  dominated  by  the  contri¬ 
bution  from  the  largest  size  interval  (6.4-12.8  pm).  For  ex¬ 
ample,  at  time  21.40  the  6.4-1 2.8-pm  channel  contributes  over 
8(Fc  of  total  extinction.  Therefore,  it  is  not  surprising  to  find 
a  strong  correlation  between  extinction  and  LWC  for  the 
6.4-12.8-pm  interval  of  Lenham  and  Clav1  since  the  channel 
dominates  both  extinction  and  liquid  water  content. 

Moreover  the  inversion  scheme2  used  by  Lenham  and  Clay1 
truncates  the  size  distribution  spectrum  at  12.8-pm  radius. 
There  is  strong  evidence8-8  that  both  radiation  and  advection 
fogs  possess  radii  >n  excess  of  12.8  pm.  The  contribution  of 
a  seventh  channel  (12.8-25.6  pm)  to  total  extinction  would 
amount  to  more  than  l(Ti  allowing  for  a  plausible  order  of 
magnitude  drop  in  concentration  cm-3  from  that  of  channel 
6.  This  wouid  grossly  alter  the  form  of  the  histogram  deduced 
by  Lenham  and  Clay,1  and  it  would  also  weaken  the  high 
correlation  between  visible  extinction  and  liquid  water  content 
obtained  by  Lenham  and  Clav.1 

The  relationship  between  calculated  exrnction  and  LWC 
based  on  measured  fog  size  distributions  for  radiation  and 
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Fig.  1.  Variation  of  extinction  coefficient  with  liquid  water  content 
in  atmospheric  fog  and  hazt  for  320  size  distribution  measurements 
made  at  dO'.erent  geographic  locales  and  under  a  variety  of  meteo- 
•ological  conditions  as  shown  by  Pinnick  et  a/.8 


advection  fog  at  X  =  0.55  /um  has  been  shown8  not  to  be 
unique.  Indeed  extinction  can  vary  by  approximately  an 
order  of  magnitude  for  a  particular  LWC  value  as  shown  by 
Pinnick  et  al.8  and  reproduced  in  Fig.  1. 

The  correlation  of  number  concentration  w-ith  extinction 
was  examined  by  Lenham  and  Clay.1  In  general,  they  found 
a  weak  correlation  between  number  concentration  and  ex¬ 


tinction  with  the  exception  of  the  6.4-12.8-/zm  category  for 
a  dense  fog.  However,  a  study  of  inferred  droplet  concen¬ 
tration  for  the  6.4-12.8-/im  category  and  extinction  at  X  =  0.53 
^m  for  the  four  fog  case  studies  of  their  Fig.  1  reveals  that  a 
strong  correlation  does  not  exist  between  number  concen¬ 
tration  and  extinction.  It  is  clear  from  Eq.  (1)  that  extinction 
will  be  linearly  related  to  number  concentration  only  if  ex¬ 
tinction  cross  section  <rext  and  size  distribution  n(r)  vary  with 
droplet  size  over  the  whole  size  range  in  an  inverse  manner. 
For  example,  if  aest  varies  as  r2,  extinction  will  vary  linearly 
with  number  concentration  only  if  n(r)  varies  as  1/r2.  The 
above  conditions  on  aext  and  n(r)  are  unlikely  to  be  met  over 
a  typical  fog  droplet  size  distribution  in  view  of  la)  the  char¬ 
acteristic  resonance  response  of  crext  to  size  and  (b)  the  wide¬ 
spread  evidence6-8  that  many  fog  droplet  size  distributions 
are  bimodal  in  character.  Thus,  in  general,  it  is  extremely 
unlikely  that  extinction  will  be  linearly  related  to  number 
concentration  for  a  fog  size  distribution  spectrum  at  visible 
wavelengths. 

This  work  was  supported  by  means  of  a  research  contract 
with  the  U.S.  Army  through  its  European  Research  Office, 
London. 

References 

1.  A.  P.  Lenham  and  M.  R.  Clay,  Appl.  Opt.  21, 4191  (1982). 

2.  J.  Heintzenberg,  H.  Mueller,  H.  Quenzel,  and  E.  Thomalla,  Appl. 
Opt.  20,1308(1981). 

3.  Petr  Chylek,  J.  Atmos.  Sci.  35,  296  (1978). 

4.  W.  J.  Wiscombe.  Appl.  Opt.  19,  1505  (1980). 

5.  G.  M.  Hale  and  M.  R.  Querry,  Appl.  Opt.  12,  555  (1973). 

6.  J.  A.  Garland.  Q.  J.  R.  Meteorol.  Soc.  97,  483  (1971). 

7.  R.  J.  Pilie,  E.  J.  Mack,  W.  C.  Kocmond,  W.  J.  Eadie.  and  C.  W. 
Rogers,  J.  Appl.  Meteorol.  14,  364  (1975). 

8.  R.  G.  Pinnick.  S.  G.  Jennings.  P.  Chylek.  and  H.  J.  Auvermann,  J. 
Atmos.  Sci.  36,  1577  (1979). 


*  September  1983  Vol  22.  No  17  /  APPLIED  OPTICS 


2515 


